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The  results  of  analytical  studies  are  presented  in  which  (l)  the  SAP  IV 
elastic  finite— element  (FE)  code  was  used  to  perform  modal  analyses  for  a 
rectangular  structure  in  vacuo  and  embedded  in  different  backfill  materials’, 
and  (2)  a series  of  dynamic  and  static  two-dimensional  (2D),  plane-strain, 

FE,  structure— medium  interaction  (SMI)  code  calculations  were  performed  using 
the  HONDO  code.  Uie  HONDO  code  SMI  calculations  were  made  to  investigate  the 
effect  of  variations  in  (l)  surface  airblast  loadings;  and  (2)  constitutive  — 
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20.  ABSTRACT  (Continued). 

^ properties  of  the  backfill  region  around  a hypothetical,  shallow-buried 

structure  on  the  transient  loading  environments  experienced  by  the  structure. 

The  modjh.  analyses  showed  that  the  natural  frequencies  increased  as 
the  quality  of  the  material  surrounding  the  structure  increases.  They 
also  indicated  that  the  planned  SMI  calculations  would  have  sufficient 
frequency  response.^ 

A first  series  of  fi>ur  dynamic  2D,  11  MI  calculations  were  performed 
on  this  pline-strain  idealization  of  a simple  buried  structure.  Thesi 
calculations  involved  tWo  surface  loading^  (o.g.,  nominal  300-psi  (2^).7-har) 
ove rpressu/re  pulses  for / 1-MT  and  10-KT’  weapons)  over  two  backfill  conditions 
(i.e.,  loose  and  dense  backfill  cases).  /A  fifth  dynamic  calculation, 
essentially  identical  to  the  first  four  calculations  withl a nominal/ 

1000-psi  (68.9-bar)  overpressure  pulse  mv  a 10-KT  weapon)  over  tlie  dense 
backfill  (condition  was  conducted.  A filial  calculation  (aJ  "pseudostjatic" 
calculation)  was  conducted  in  which  the  structure  surrounded  by  dcrjse 
backfill; was  loaded  with  a monotonicaliy  increasing  load /from  0 to j300  psi 
(0  to  20 j. 7 bars). 

/ ! 

Comparisons  of  insults  of  the  firdt  series  of  ealcu  ations  indicated 
that  thb  use  of  both  loose  backfill  rather  than  dense  backfill  and  higher 
weapon-jyield  loading  oases  increased  (1)  deflections  across  backfill 
sections;  (2)  loads  on,  deflections  of,  and  thrusts,  shears,  and  bending 
momenta  within  the  structure;  and  (3)  Amplitudes  of  the  shock  spectra  for 
points  on  the  inside  Surface  of  the  structure.  The  results  of  the  fifth 
calculation  indicated  that  the  stressed  and  bending  momfnt  within  and  the 
deflections  of  the  stfucture  for  the  nominal  1000-psi  (68.9-bar)  overpressure 
airblast  for  a 10-KT  Weapon  over  the  dense  backfill  problem  were  almost 
identical  to  those  for  the  300-psi  (20.7-bar)  airblast  (for  a 1-MT  weapon 
over  the  dense  backfill  problem. 

V There  was  a noticeable  reduction  in  the  frequencies  at  which  the 
structures  responded  in  the  HONDO  code  calculations  due  to  a decrease  in 
stiffness  of  the  soil  backfill.  A linear  elastic  single-degree—.' f- freedom 
system  model  based  on  the  predominant  natural,  mode  of  the  embedded  structure 
simulated  the  response  of  the  structure  to  the  three  different  airblast 
loadings  only  when  a 20  to  to  percent  of  critLcal  damping  was  applied.  The 
damping  is  believed  to  be  necessary  to  account  for  energy  removed  from  the 
elastic  structure  because  of  radiation  into  the  surrounding  media  and  the 
hysteretic  nature  of  this  media. 
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EFFECT  OF  BACKFILL  PROPERTY  AND  AIRBLAST  VARIATIONS  ON 
THE  EXTERNAL  LOADS  DELIVERED  TO  BURIED  BOX  STRUCTURES 


CHAPTER  1 
INTRODUCTION 


1.1  BACKGROUND 

Large  amounts  of  backfill  must  be  placed  around  and  over  buried 
hardened  facilities  for  hardening  against  nuclear  airblast  and  ground 
shock  environments.  The  condition  of  the  backfill  around  present  fa- 
cilities ranges  from  very  veil  compacted  to  poorly  compacted  backfill . 
Reference  1 presents  the  results  of  a series  of  two-dimensional  (2D), 
dynamic  finite  element  (FE),  structure-medium  interaction  (SMI)  code 
calculations , which  were  conducted  to  parametrically  investigate  the 
effect  of  variations  in  constitutive  properties  of  the  backfill  region 
around  a hypothetical  plane-strain  idealization  of  a simple,  shallow- 
buried,  box-type  protective  structure.  In  this  reference,  the  differ- 
ences in  the  dynamic  response  of  the  structure  due  to  a long-duration, 
50-psi  ( 3.1*5-bar ) surface  airblast  loading*  caused  by  changing  the  sur- 
rounding backfill  from  a dense  (or  well-compacted)  glacial  till  to  the 
same  material  in  a loose  (or  poorly  compacted)  condition  were  examined. 

A second  study,  identical  with  the  first  except  that  dense  and  loose 
clay  shale  materials  were  simulated,  is  also  reported. 

Hie  calculations  performed  in  Reference  1 indicated  that  for  the 
particular  idealized  problem  investigated,  the  use  of  loose  backfill 
rather  than  dense  backfill  results  in  (1)  increased  vertical  deflections 
in  the  backfill,  (2)  increased  loads  on,  deflections  of,  and  thrusts, 
shears,  and  bending  moments  within  the  structure,  and  (3)  increased 
amplitudes  of  the  shock  spectra  for  points  on  the  inside  surface  of  the 
structure . 

Because  of  the  long  airblast  positive-phase  duration,  negligible 


* Due  to  a 1-MT  surface  burst. 
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stress  attenuation  occurred  in  the  loose  backfill  calculation  in  spite 
of  the  high  hysteresis  of  that  material.  This  should  not  be  the  case 
for  short-duration  airblast  Loadings. 

Under  the  present  study.  Defense  Nuclear  Agency  (DNA)  Subtask  SC062, 
Work  Unit  01,  "influence  of  Backfill  on  Structural  Response,"  five  2D 

dynamic  and  one  2D  static  FE,  SMI  code  calculations  were  performed  with 

2 i 

the  HONDO  * code  to  determine  the  effect  of  soil  backfill  property 
variations  and  surface  airblast  loading  options  on  the  transient  ground 
shock  environments  experienced  by  a hypothetical  protective  structure 
at  much  higher  overpressure  levels  than  those  studied  in  Reference  1. 

The  results  are  primarily  intended  to  (1)  determine  the  effect  of  soil 
property  variations  and  surface  loading  options  on  the  transient  ground 
shock  environments  experienced  by  a hypothetical  shallow-buried  structure, 
and  (2)  provide  a tangible  example  to  structure  analysts  of  the  type  of 
structure/backfill  interaction  information  that  can  be  obtained  for 
their  use  in  vulnerability  assessments. 

Prior  to  performing  these  calculations,  a series  of  in  vacuo  modal 

4 

analyses  was  conducted  with  the  SAP  IV  elastic  PE  code  to  (1)  determine 
the  effects  of  using  basic  isoparametric  quadrilateral  elements  (i.e., 
the  type  employed  by  the  HONDO  code),  instead  of  high  order  elements, 
on  calculated  natural  frequencies,  and  (2)  aid  in  determining  the  mini- 
mum number  of  elements  needed  across  structural  sections  in  order  to 
adequately  model  the  response  of  the  structure,  Embedded  modal  analyses 
were  also  performed  to  (1)  determine  the  effect  of  material  property 
variations  on  the  natural  frequencies  of  the  buried  structure,  and  (2) 
assist  in  planning  the  follow-on  structure/backfill  interaction  calcula- 
tions. The  results  of  the  modal  analyses  were  previously  submitted  in 
a letter  report  to  DNA*  and  are  included  herein  as  Appendix  A. 


J.  0.  Curtis  and  J.  E.  Windham,  "in  Vacuo  and  Embedded  Modal  Analyses 
for  Hypothetical  Structure  in  Plane  Strain,"  SC002/01  informal  report 
to  HQ  DNA,  April  1976,  U.  S.  Army  Engineer  Waterways  Experiment  Station, 
CE,  Vicksburg,  Mississippi. 
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A first  series  of  HONDO  code  calculations  consisting  of  four  large 
2D,  SMI  calculations  were  performed  on  a GDC  7600  computer.  These  calcu- 
lations involve  tvo  surface  loadings  (e.g.,  nominal  300-psi  (20.7-bar) 
overpressure  pulses  for  1-MT  and  10-KT  weapons)  over  two  backfill  condi- 
tions (i.e.  loose  and  dense  shale  backfill  cases'1').  A detailed  plan  of 
study  for  this  first  series  of  calculations  has  previously  been  submitted.* 
A fifth  dynamic  calculation,  essentially  identical  to  the  other  dense 
shade  backfill  cadculations  except  for  the  surface  blast  loading  that 
was  a nominal  1000-psi  (68.9-bar)  overpressure  pulse  from  a 10-KT  detona- 
tion, was  also  conducted.  Finedly,  a calculation  was  conducted  with  the 
HONDO  code  in  which  the  structure  surrounded  by  dense  backfill  was  slowly 
loaded  with  a mono tonic ally  increasing  load  from  0 to  300  psi  (20.7  bars). 
The  stress  rate  was  applied  slow  enough  so  that  the  results  at  any  stress 
level  would  closely  approximate  that  of  a static  solution.  In  order  to 
maximize  correlation  with  previous  results,  the  site  profile,  excavation, 
and  structure  configurations  were  identical  to  those  in  Reference  1. 

1.2  PURPOSE  AND  SCOPE 

The  purposes  of  this  code  calculation  study  are  to  (l)  determine  the 
effect  of  selected  soil  backfill  property  variations  and  surface  loading 
options  on  the  high  overpressure,  airblast- induced  ground  shock  environ- 
ments experienced  by  a hypothetical  shallow-buried  structure,  and  (2) 
provide  a tangible  example  to  structural  analysts  of  the  type  of  backfill- 
structure  interaction  information  that  can  be  obtained  for  use  in  buried 
structure  vulnerability  assessments.  This  report  (l)  outlines  the  plan  of 
analyses,  (2)  describes  the  FE  code  calculations,  (3)  presents  their 
results,  (U)  reports  comparative  analyses  of  the  effects  of  the  indepen- 
dent variables,  and  (5)  discusses  the  implications  of  these  results. 

The  analyses  are  essentially  limited  to  the  results  of  the  six  HONDO  code 


* J.  E.  Windham,  "Plan  of  Study  for  First  Series  of  Structure /Backfill 
Interaction  Calculations,"  WES  informal  report  on  DNA  Subtask  SC0O2, 
Work  Unit  01,  April  1976. 


| 


1 


U 


15 


calculation:)  that  are  the  key  elements  of  thi3  study. 

A supplementary  study  is  reported  in  Appendix  B.  A simple  method 
for  producing  the  rigid  body  motion  of  a buried  box  structure  is  pre- 
sented in  Keference  5.  This  method  is  applied  to  the  1-MT,  300-psi 
(20. 7-bar)  dense  backfill  case.  The  calculated  vertical  motion  of  the 
structure  was  found  to  be  in  good  agreement  with  the  HONDO  code  calcu- 
lated structure  motion  time  histories. 
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CHAPTER  2 


DESCRIPTION  OF  FIRST  CALCULATION  SERIES 


2.1  OVERALL  PLAN 


The  objective  of  the  first  series  of  calculations  was  a parametric 
study  to  determine  the  effects  of  both  (1)  soil  backfill  property  vari- 
ations (due  to  either  diverse  construction  conditions  or  data  uncer- 
tainties), and  (2)  surface  airblast  loading  options  (due  to  selection 
of  different  weapon  yields)  on  the  transient  ground  shock  environments 
experienced  by  a hypothetical  protective  structure  located  at  a given 
overpressure  contour.  In  order  to  maximize  correlation  with  previous 
results,  the  site/excavation/structure  configuration  is  identical  with 

that  reported  in  Reference  1.  A total  of  four  large  calculations  were 

2 3 

performed  on  a CDC  7600  computer  with  the  HONDO  * dynamic  FE  code; 
these  calculations  involve  two  surface  loadings  (e.g.,  nominal  300-psi 
(20.7-bar)  overpressure  pulse  for  1-MT  and  10-KT  weapons)  over  each  of 
two  different  backfills  (loose  and  dense  shale). 

Figure  2.1  illustrates  the  geometry  of  the  plane-strain  problem  in- 
vestigated. The  free-field  medium  is  postulated  ac  a four-layered  soil 
profile  consisting  of  clayey  sand  over  clay  shade  (Faterial  No.  1 through 
4).  Ihe  structure  is  assumed  to  be  constructed  of  reinforced  concrete 
(Material  No.  5)*  which  is  idealized  as  a linear  elastic  material.  It  is 
laterally  supported  and  covered  by  a compacted  soil  backfill  (Material 
No.  6),  whose  constitutive  properties  are  independent  variables  in  this 
study. 


2.2  FINITE  ELEMENT  GRID 


'JBie  FE  grid  that  was  employed  is  shown  in  Figure  2.1  and  consisted 
of  1048  constant  strain  rectangular  and  triangular  elements.  Hie 
elements  located  in  the  structure  are  1.0  foot  (0.30  metre)  by  1.25  feet 
(0.38  metre),  1.0  foot  (0.30  metre)  by  2.5  feet  (0.?6  metre),  and 
1.25  feet  (0.38  metre)  by  2.0  feet  (0.6l  metre);  each  structural  member 


has  four  elements  across  its  section.  In  the  backfill,  the  elements 
range  from  1.25  feet  (0.  jS  metre)  by  2.5  feet  10. '(6  metre)  to  2.5  feet 
(0.76  metre)  by  10  feet  (3.05  metres).  In  the  stiffer  and  deeper  earth 
materials,  the  elements  are  larger.  'lhe  largest  elements  are  in  the 
corners  of  the  grid  and  are  20  feet  (c.10  metros)  by  20  feet  (6.10  metres). 
The  grid  was  chosen  as  a compromise  between  (l)  the  desire  for  fine  reso- 
lution, and  (2)  the  need  to  keep  computer  memory  and  running  time  (and, 
hence,  cost)  requirements  within  reason. 

2.3  HYPOTHETICAL  STRUCTURE 

The  idealised  structure  is  covered  by  5 feet  (1.52  metres)  of  back- 
fill and  is  supported  on  undisturbed  shale  at  a depth  of  30  feet 
(9.1U  metres);  th  structure  has  outside  dimensions  of  **0  feet 
(12.19  metres)  and  25  feet  (7.b2  metres).  Its  roof  and  floor  are 
5 feet  (1.52  metres)  thick.  Its  sidewalls  are  U feet  (1.22  metres)  thick. 
The  exterior  structure  elements  ore  treated  as  being  fully  bonded  to 
those  of  the  backfill  and  underlying  shale,  as  no  slip  element  capability 
existed  in  the  U.  S,  Army  Engineer  Waterways  Experiment  Station  ( WES ) 
version  of  the  HONDO  code  when  these  calculations  were  performed.*  The 
structure  is  treated  as  a linear  elastic  material  with  a bulk  modulus  K 
of  1330  ksi  (91,700  bars)  and  a shear  modulus  0 of  SOO  ksi  (55,200  bars ) . 
The  structure  is  assumed  to  have  a mass  density  of  1^*5  pcf  (2.12  g/em^). 
Whether  or  not  a real  reinforced  concrete  structure  of  the  dimension 
adopted  would  remain  operational  in  the  overpressure  environments  used 
in  this  study  is  doubt  Ail.  However,  since  (l)  the  purpose  of  the  planned 
calculations  was  to  show  only  the  relative  effect  of  weapon  yield  on  the 
loads  transmitted  to  a liypothetlcal  structiu'e  for  a wide  range  of  backfill 
materials,  and  (2)  it  was  desirable  for  comparison  purposes  to  use  the 
same  structure  as  that  used  in  previous  calculations,1  the  decision  was 
made  to  use  the  section  thicknesses  indicated  in  Figure  2.1. 

* A slip  de bonding  line  capability  has  recently  been  added  to  the  code. 
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2.U  CONSTITUTIVE  PROPERTIES 

All  of  the  earth  materials  were  represented  in  the  calculations  by 
nonlinear  hysteretic  soil  cap  models  of  the  type  used  by  the  Site  Defense 
Ground  Motion  Working  Group  ( SDGMWG ) . 0 However,  for  the  SDGMWG  analyses, 
uniaxial  strain  (UX)  relations  were  only  modeled  to  a vertical  stress  o 

£ 

of  1000  psi  (68.9  bars).  For  the  current  study,  it  was  necessary  to 

extend  the  cap  model  fits  to  a o in  UX  of  5000  psi  (3M.7  bars).  Ihe 

z 

extended  constitutive  model  fits  for  the  in  situ  or  free-field  materials 
(Material  Nos.  1 through  4)  are  shown  in  Figures  2.2  through  2.9;  values 
of  the  fitting  constants  that  were  actually  used  as  calculational  input 
are  presented  in  Ihble  2.1. 

Ihe  backfill  material  properties  chosen  for  this  study  were  the  same 

as  those  previously  used  to  represent  a crushed  shale  backfill  placed  in 

very  loose  and  very  dense  conditions.  'Ihe  cap  model  fits  for  the  loose 

and  dense  backfill  materials,  which  are  extended  from  160  psi  (11. 0 bars) 

to  5000  psi  (3M.7  bars)  for  the  current  study,  are  shown  in  Figures  2.10 

through  2.13;  values  of  the  fitting  constants  are  also  listed  in  Table  2.1. 

The  UX  0 versus  vertical  strain  o relations  for  the  loose  and  dense 
z z 

backfills  are  compared  in  Figure  2.lU.  The  initial  constrained  modulus 
for  the  loose  backfill  material  is  1500  psi  (103.1*  bars)  and  for 
the  dense  material  is  6000  psi  (413.7  bars).  The  secant  moduli  to 
300  psi  (20.7  bars)  for  the  loose  and  dense  backfill  materials  are 
2500  psi  (172.U  bars)  and  7500  psi  (517.2  bars),  respectively. 

A review  of  UX  data  from  tests  conducted  on  backfill  materials  from 
19  sites  (including  M1NUTEMAN,  HEST,  and  high-explosive  test  event  sites) 
showed  that  the  average  secant  loading  modulus  to  300  psi  (20.7  bars) 
ranged  from  5000  psi  (3M.7  bars)  to  13,900  psi  (958.6  bars).  Backfill 
density  control  specifications  were  in  effect  at  all  of  these  sites.  Thus, 
the  dense  shale  properties  proposed  for  use  in  this  study  are  clearly  with- 
in the  range  produced  by  good  construction  practice.  However,  it  is  also 
clear  that  the  dense  backfill  should  be  considered  as  ’’typical"  rather  than 
as  an  upper  bound  to  the  backfill  stiffness  that  can  be  obtained.  Hie 
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loose  shale  properties  are  considered  to  be  typical  of  poor  construction 
practice  or  expedient  construction. 

2.5  SURFACE  OVERPRESSURE 

iVo  traveling  surface  airblast  overpressure  functions  were  used  in 

the  calculations.  The  loading  portions  for  both  are  defined  by  a linear 

rise  to  peak,  pressure  in  a constant  time  of  10  ms;  the  decay  portions 
7 

are  fits  to  Erode' s equations  for  10-KT  and  1-MT  weapons  detonated  at  a 
sero  height  of  burst  and  at  a distance  from  the  structure  that  will 
cause  a 300-psi  (20.7-bar)  surface  overpressure  over  the  centerline  of 
the  structure.  The  variations  of  the  pressure  time  histories  for  both 
cases  as  the  airblast  traverses  the  2*»0-foot  (73.2-metre)  distance  from 
the  left  to  the  right  boundary  of  the  problem  are  shown  in  Figure  2.15. 

The  pressure  time  histories  over*  the  left  boundary,  the  structure  center- 
line,  and  right  boundary  for  the  two  weapon  yields  are  plotted  at  larger 
scales  for  comparison  purposes  in  Figures  2.16,  2.17 » and  2.18,  respec- 
tively. Artificial  10-ms  rise  times  appended  to  the  Brode  overpressure 
pulses  increase  the  impulse  at  the  300-psi  (20.7-bar)  level  by  2.8  percent 
for  the  1-MT  case  and  12.9  percent  for  the  10-KT  case;  at  higher  over- 
pressure levels,  the  added  percentage  Is  even  larger.  Based  on  criteria 
developed  for  the  limiting  frequencies  that  can  be  propagated  in  ID 
elastic  wave  propagation  problems,0  even  longer  rise  times  should  be 
employed  if  the  lowest  moduli  for  the  backfill  materials  are  used  to 
calculate  the  elastic  wave  speeds  and  the  FE  grid  remains  as  defined  in 
Figure  2.1.  However,  appending  a longer  rise  time  to  the  airblast 
functions  would  unrealistically  distort  the  Impulse.  The  desired  frequency 
resolution  could  be  fully  satisfied  by  increasing  the  number  of  elements 
used  to  sone  the  problem;  but  this  would  be  cost-prohibitive.  Thus,  a 
10-ms  rise  time  is  a compromise  value,  which  may  be  artificially  increased 
by  the  FE  grid  after  the  stress  wave  has  propagated  a short  distance 
into  the  backfill  media. 

The  duration  of  each  of  the  calculations  was  200  ms . At  200  ms , the 


P 


surface  overpressures  for  the  10-KT  and.  1-MI  loadings  were  approximately 
5 psi  (0.31*  bars)  and  65  psi  (4.48  bars),  respectively.  Reflections 
from  the  bottom  and  right  boundaries  reached  the  structure  at  a time  of 
approximately  100  ms.  After  that  time,  there  was  some  influence  of  the 
boundary  on  structure/backfill  interaction,  but  the  influence  should  be 
nearly  the  same  in  all  four  problems  since  the  first  reflected  signals 
travel  primarily  through  the  nonvaried  free-field  materials. 

2.6  OUTPUT  DATA 

Hie  output  data  from  the  calculations  consisted  of  (l)  time  histories 
of  stress  and  motion  at  selected  locations  in  the  structure  and  in  the 
earth  media,  (2)  plots  of  the  normal  stress  on  the  exterior  surfaces  of 
the  structure,  (3)  plots  of  deflected  shape  of  the  structure's  neutral 
axis,  and  (4)  tabulated  bending  moment,  axial  thrust,  and  shear  at 
selected  structure  cross  sections  for  selected  times.  Hie  locations  for 
which  time  histories  of  motion  and  stress  were  saved  are  shown  in  Fig- 
ure 2.19.  Approximately  260  motion  and  200  stress  time  histories  were 
involved  for  each  calculation.  Forty-seven  plots  (or  snapshots)  of 
normal  stress  on  the  external  surface  of  the  structure  and  of  the  de- 
flected shape  of  the  structure,  and  the  tables  of  bending  moment,  axial 
thrust,  and  shear  were  made  at  2. 5-ms  time  increments  between  10  and 
100  ms  and  at  10-ms  increments  between  100  and  200  ms.  The  28  structural 
sections  for  which  bending  moments,  axial  thrusts,  and  shear  were  tabu- 
lated are  indicated  in  Figure  2.20. 


2.7  TIME  INCREMENT  AND 

FREQUENCY  TRANSMISSION 

Each  calculation  was  run  for  2000  time  steps  of  0.1  ms  each.  This 
time  step  was  chosen  to  satisfy  the  Courant  criteria*  and  was  controlled 


. ..  . min  . 1.0  ft  . 

i — i 8.76  rt/«s  i °-Ul‘  “ 


where  At  = time  increment;  Ax 


minimum  element  size;  Cp 


maximum  P-wave  velocity. 
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by  the  minimum  FE  dimension  and  P-wave  velocity  within  the  elastic  struc- 
ture. The  lowest  frequencies  that  can  be  fully  transmitted  are  given  by 
the  reciprocal  of  the  calculation  pulse  duration,  i.e.,  1/200  ms  or  5 Hz. 
The  lowest  modal  frequency  determined  from  an  embedded  modal  analysis 
(see  Appendix  A)  using  lower  bound  moduli  for  all  materials  was  4.31$  Hz, 
while  a similar  analysis  using  upper  bound  moduli  gave  a first-mode  fre- 
quency of  11.99  Hr.  Thus,  the  duration  of  the  calculations  were  long 
enough  to  permit  peak  displacement  responses  in  the  lowest  mode  of  the 
system.  The  highest  frequencies  that  can  be  fully  transmitted  cannot 
exceed  the  reciprocal  of  the  loading  rise  time,  i.e.,  1/10  ms  or  100  Hz. 
As  a practical  matter,  however,  credible  frequencies  probably  did  not 
exceed  one-half  of  this  value  or  50  Hz. 

The  FE  grid  for  the  backfill  material  immediately  surrounding  the 
structure  can  transmit  frequencies  through  the  dense  backfill  up  to 
68  Hz  based  on  the  formula 


max  wAx 


(2.1) 


By  similar  calculation,  however,  the  loose  backfill  can  only  transmit  fre- 
quencies up  to  24  Hz.  Thus,  the  cutoff  frequencies  were  assumed  to  be 
on  the  order  of  50  Hz  for  the  two  dense  backfill  calculations  and  25  Hz 
for  the  two  loose  backfill  calculations.  The  embedded  modal  analyses 
presented  in  Appendix  A using  upper  bound  moduli  gave  a 20th-raode  fre- 
quency of  49,11  Hz;  the  calculated  frequency  for  the  25th  mode  based  on 
lower  bound  moduli  was  only  14.35  Hz.  Thus,  the  calculations  were  Judged 
to  have  a sufficient  range  of  frequency  response  to  permit  an  adequate 
assessment  of  the  stress  and  motion  patterns  that  would  occur  within  the 
hypothetical  structure/backfill  system  under  the  postulated  airblast 
environments. 
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try  of  problem  and  FE  grid  representation 


PRINCIPAL  STRESS  DIFFERENCE 
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Figure  2.3  Failure  envelope  and  UX  stress  path  for  free-field 
layer  1 produced  by  the  constitutive  model. 
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Figure  2.5  Failure  envelope  and  UX  atreea  path  for  free-field  layer  2 
produced  by  the  constitutive  model. 
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TIME  AFTER  BEGINNING  OF  CALCULATION  (I  MT  LOADING),  MSEC 


Figure  2.18  Surface  overpressure  time  history  over  the  right  boundary 
of  the  problem  for  1-MT  and  10-KT  weapons. 
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Figure  2.19  Detailed  layout  of  the  FE  grid  in  the  vicinity 
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RESULTS  FROM  AND  ANALYSIS  OF  1HE 
FIRST  SERIES  OF  COMPUTATIONS 

3.1  STRESS  AND  DISPLACEMENT  PATTERNS 

AT  SELECTED  TIMES 

3.1.1  Deflection  Across  Backfill  Sections 

Calculated  vertical  displacement  patterns  for  sections  through  the 
loose  and  dense  backfills  at  60  and  100  ms  into  the  calculations  sure  shown 
in  Figures  3.1  and  3.2,  respectively,  for  the  1-MT  loading  case  and  in 
Figures  3.3  and  3.4  for  the  10-KT  loading  case.  As  can  be  seen  in  all 
cases,  the  vertical  deflections  in  the  loose  backfill  at  the  5-foot 
(1.52-metre)  depth  are  much  larger  at  60  ms  than  those  in  the  dense 
backfill.  This  appears  reasonable  when  the  stress-strain  curves  for 
the  two  materials  are  examined  (Figure  2.14).  But  displacements  at  the 
15-foot  (4.57  metre)  depth  in  the  dense  backfill  at  this  time  axe  larger 
than  those  for  the  loose  backfill  because  the  higher  wave  velocity  of 
the  dense  shale  fill  has  allowed  the  peak  stress  to  propagate  to  a 
greater  depth.  At  100  ms,  however,  when  wave  propagation  no  longer 
plays  this  type  of  interactive  role  and  the  peak  free-field  stress  has 
occurred  at  all  depths  of  interest,  the  displacements  are  everywhere 
larger  in  the  loose  backfill  than  in  the  dense  backfill,  as  shown  in 
Figures  3.2  and  3.4,  respectively.  At  60  ms,  the  vertical  backfill 
deflections  for  the  1-MT  loading  case  are  only  slightly  higher  than 
those  for  the  10-KT  loading;  at  100  ms,  the  vertical  deflections  are 
much  higher  for  the  1-MT  loading. 

The  maximum  deflections  of  the  backfill  (node  148)  at  the  5- foot 
( 1.52-metre)  depth  and  the  times  that  they  occur  along  with  the  strain 
at  200  psi  (20.7  bars)  from  the  UX  stress-strain  curves  for  the  two 
backfill  materials  are  shown  below: 
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Deflection  at  5- ft  (1.52-m)  Depth 

Node  148 

1-MT  Weapon  Yield  10-KT  Weapon  Yield 


Peak  Deflection  Time  Peak  Deflection  Time 
cm ms  cm  ms 


Dense 

Loose 


3.6 
13.  4 


24  85  18  88 

76  120  53  115 

As  can  be  seen  in  the  table,  the  deflection  at  the  5-foot  (1. 52-metre) 
depth  increases  as  weapon  yield  increases  and  as  backfill  quality  de- 
creases. For  each  weapon  yield,  the  deflection  at  the  5- foot  (1.52-metre) 
depth  for  the  loose  backfill  case  is  approximately  three  times  that  for  the 
dense  backfill  case.  The  deflections  in  the  backfill  for  the  1-MT  weapon 
yields  over  loose  and  dense  backfills  are  a factor  of  1.4  and  1.3  times 
the  deflections  for  the  10-KT  weapon  yields  over  loose  and  dense  backfills, 
respectively.  These  trends  are  reasonable  when  the  stress-strain  proper- 
ties of  the  backfill  materials  and  the  differences  in  propagation  speeds 
and  impulse  for  the  two  airblast  loadings  are  considered. 

'Hie  interface  between  the  backfill  and  structure  should  have  been 
treated  computationally  by  some  type  of  slip  element  so  that  the  rela- 
tively compressible  backfill  could  deflect  vertically  with  respect  to 
the  concrete  structure  (which  is  founded  on  very  stiff  undisturbed  clay 
shale).  This  type  of  element  was  not  available  when  these  calculations 
were  performed.  Since  the  structure  and  backfill  were  "welded"  together, 
it  is  inevitable  that  very  little  deflection  would  occur  near  the  struc- 
ture walls.  Of  course,  this  is  not  the  case  in  the  real  world.  However, 
the  backfill  vertical  deflections  did  maximize  at  a distance  of  approxi- 
mately 6 to  10  metres  from  the  structure  (Figures  3.1  through  3.4)  for 
all  cases.  Deflections  close  to  but  perhaps  less  than  the  maximum  shown 
for  the  backfill  (because  of  wall  friction)  should  be  expected  close  to 
the  structure,  lhe  deflections  calculated  in  the  loose  shale  backfill 
would  be  considered  very  severe  for  cables  or  pipes  that  might  be  con- 
nected to  the  structure  at  some  point  along  the  sidewall. 
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3.1.2  Stresses  on  the  Exterior  Structure  Siarfaces 
and  Deflections  of  the  Neutral  Axis 

Figures  3.5 » 3.7»  3.8,  3.9,  and  3.10  show  the  stresses  normal  to  the 
exterior  surfaces  of  the  structure  and  the  deflection  of  the  neutral  axis 
of  the  structural  members  for  all  four  calculations  at  times  of  U0 , 80, 

100,  150,  and  200  ms,  respectively.  Figure  3.6  shows  the  same  normal 
stress  and  deflection  information  for  the  times  that  the  vertical  deflection 
at  the  center  of  the  roof  of  the  structure  has  maximized  for  each  calcula- 
tion. Maximum  deflections  occur  between  60  and  70  ms  for  all  calculations. 
In  actuality,  the  stresses  are  the  vertical  or  horizontal  stresses  in  the 
middle  of  the  soil  "cells"  adjacent  to  the  structure. 

At  a time  of  40  ms  as  shown  in  Figure  3.5,  the  roof  loading  is  tri- 
angular except  in  the  10-KT  dense  backfill  case.  At  40  ms,  the  peak  over- 
pressure for  the  10-KT  dense  backfill  simulation  has  already  been  applied 
to  the  structure  roof,  while  the  stress  in  the  other  cases  is  still  in- 
creasing. The  amplitude  of  the  stress  normal  to  the  roof  at  this  time  is 
a function  of  the  velocity  of  the  traveling  airblast  and  the  wave  speed 
of  the  backfill.  The  average  propagation  velocity  of  the  10-KT  airblast 
for  locations  upstream  of  the  center  of  the  structure  is  faster  than  that 
of  the  1-MT  airblast  (i.e.,  the  local  propagation  velocity  is  2.0  m/ms 
versus  1.51*  m/ms  at  the  left  (upstream)  boundary  of  the  problem),  and  the 
wave  speed  of  dense  backfill  (0.195  m/ms)  is  greater  that  that  of  the 
loose  backfill  (0.115  m/ms).  At  the  time  that  the  peak  centerline  deflec- 
tions of  the  roof  occur,  the  roof  loading  varies  irregularly  with  position, 
but  generally  the  differences  among  the  four  calculations  are  less  them 
20  bars  (Figure  3.6).  For  subsequent  times,  i.e.,  80,  100,  150,  and 
200  ms  as  shown  in  Figures  3.7  through  3.10,  the  normal  stress  on  the 
roof  decreases  with  increasing  time  and  tends  toward  a pattern  of  uniform 
loading  over  most  of  the  roof  with  stress  concentrations  a^ove  the  side- 
walls.  At  200  ms,  the  stresses  normal  to  the  roof  for  the  1-MT  weapon- 
yield  cases  are  higher  than  those  for  the  10-KT  cases  because  the  airblast 
from  the  10-KT  weapon  decays  at  a faster  rate  than  does  the  airblast  from 
the  1-MT  weapon. 
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The  normal  stresses  on  the  blastward  sidewall  from  the  two  dense 
backfill  calculations  at  40  msec,  as  shown  in  Figure  3.5,  are  higher 
them  those  for  the  loose  backfill  calculations;  the  higher  wave  velocity 
of  the  dense  backfill  materials  has  permitted  significant  stress  to  propa- 
gate further  down  the  backfill.  Figure  3.6  also  shows  this  trend.  At 
80  ms,  as  shown  in  Figure  3. 7.  the  stress  normal  to  the  blastward  side- 
wall  is  nearly  equal  for  both  backfill  cases. 

At  U0  ms,  the  stress  normal  to  the  structure  floor  shows  a concentra- 
tion under  the  blastward  sidewall  for  all  cases  because  significant  stress 
has  just  propagated  to  this  depth  on  the  blastward  side  of  the  structure. 
The  load  distribution  shown  in  Figure  3.5  causes  a rotation  in  the 
counter-clockwise  direction.  At  the  time  that  the  peak  deflection  of 
the  roof  is  occurring,  the  stress  normal  to  the  floor  of  the  structure 
is  concentrated  under  the  sidewalls  and  is  at  a minimum  near  the  center. 

For  subsequent  times,  as  shown  in  Figures  3.7  through  3.10,  the  stresses 
under  the  sidewalls  decrease.  At  150  ms,  as  shown  in  Figure  3-9,  tensile 
stresses  are  occurring  under  the  left  portion  of  the  floor  of  the  struc- 
ture. It  is  separating  from  the  soil,  and  a clockwise  rotation  of  the 
structure  is  taking  place  at  this  time. 

The  deflected  shape  of  the  structure  at  40  ms , as  shown  in  Figure  3.5, 
indicates  that  the  left  portion  of  the  roof  and  floor  of  the  structure  are 
deflecting  downward.  The  deflections  are  greatest  for  the  10-KT  dense 
backfill  calculation  and  the  least  for  the  1-MT  loose  backfill  calculation. 
This  is  consistent  with  the  stresses  normal  to  the  structure  at  this  time 
(Figure  3.5).  Figure  3.6  shows  the  deflected  shape  of  the  structure  at 
the  time  that  the  peak  downward  deflection  of  the  roof  centerline  occurs 
in  each  calculation.  Also,  in  Figure  3.6,  the  roof  centerline  is  deflect- 
ing downward  with  respect  to  the  sidewalls  in  all  cases,  and  the  magnitude 
of  the  deflection  increases  with  increase  in  weapon  yield  and  with  loose 
backfill.  The  total  maximum  downward  motion  of  the  roof  centerline  for 
the  1-MT  dense,  1-MT  loose,  10-KT  dense,  and  10-KT  loose  calculations  are 
10.9,  13.0,  7*2,  and  8.4  centimetres,  respectively.  The  maximum  relative 
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downward  displacements  of  the  roof  centerline  (i.e.,  measured  vertically 
from  a line  drawn  through  the  two  upper  corners  of  the  structure),  are 
5.3,  7.4,  3.4,  and  4.6  centimetres,  respectively.  The  maximum  relative 
deflections  for  the  1-MT  case  are  56  percent  and  6l  percent  greater  than 
those  for  the  10-KT  loading  cases  for  the  dense  and  loose  backfills , re- 
spectively. The  loose  backfill  caused  maximum  relative  deflections 
40  percent  higher  for  the  1-MT  loading  case  and  35  percent  higher  for 
the  10-KT  loading  case  than  did  the  dense  backfill. 

At  SO  ms,  the  roof  structure  is  rebounding  upward  with  respect  to 
the  top  corners,  while  the  structure  as  a whole  is  moving  downward 
(Figure  3.7).  The  roof  he.  rebounded  more  for  the  10-KT  loading  cases 
than  for  the  1-MT  loading  cases;  this  phenomenon  is  compatible  with  the 
stress  distributions  normal  to  the  structure  roof  shown  in  Figure  3.7. 

At  100  ms,  as  shown  in  Figure  3.8,  the  roof  has  rebounded  in  all 
cases,  and  the  downward  deflection  of  the  centerline  of  the  floor  is 
"catching  up"  with  the  sidewalls.  At  this  time,  the  structure  has 
moved  downstream  and  an  inflection  has  developed  in  the  blastward  side- 
wall  (i.e.,  the  top  of  the  wall  has  moved  further  downstream  than  the 
bottom).  At  200  ms,  as  shown  in  Figure  3.10,  for  the  10-KT  loading 
cases,  the  structure  has  essentially  returned  to  an  unde formed  condition 
but  has  been  moved  downstream  and  rotated  in  the  clockwise  direction. 

At  the  same  time,  the  roof  and  the  floor  of  the  structure  have  deflected 
inward,  and  the  sidewalls  have  deflected  outward  for  the  1-MT  loading 
cases.  The  structure  has  been  deflected  downward  further  into  the  clay 
shale  and  undergone  less  clockwise  rotation  for  the  1-MT  cases  than  for 
the  10-KT  loading  cases. 

3.1.3  Rotation  of  the  Structure 

The  rotations  of  the  chord  formed  by  the  two  bottom  outside  comers 
of  the  structure  are  presented  in  Table  3.1  for  all  calculations.  In  all 
cases,  the  structure  initially  undergoes  a counterclockwise  rotation;  then, 
as  the  airblast  passes,  it  rotates  through  the  upright  or  "rero  position" 
and  undergoes  a clockwise  rotation.  For  both  loading  conditions,  at 
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200  ms,  the  structure  remains  rotated  in  a clockwise  position.  However, 
for  the  1-MT  cases  at  200  ms,  the  clockwise  rotation  is  less  than  the 
maximum;  whereas,  for  the  10-KT  calculations  the  maximum  clockwise  rota- 
tion occurs  at  200  ms.  The  maximum  counterclockwise  rotations  are 
higher  for  the  higher  weapon  yields  and  for  loose  backfill  conditions. 

The  maximum  rotation  of  0.002U6  radians  for  the  10-KT  simulation  repre- 
sents only  3 centimetres  of  differential  displacement  over  the  length  of 
the  structure. 

3.1.U  Thrusts,  Shear  and  Bending 

Moments  Within  the  Structure 

Axial  thrust,  bending  moment,  and  snear  per  unit  width  of  the  struc- 
ture at  selected  times  were  calculated  from  the  output  data  and  were  tabu- 
lated in  Tables  3.2  through  3.10  for  the  structural  section  shown  in 
Figure  2.20.  The  sign  convention  for  these  tabulated  values  is  designated 
in  the  figure.  Shear  and  moment  diagrams  for  each  calculation  were 
plotted  for  the  roof  and  floor  at  the  times  that  the  horizontal  stress 
in  elements  511*  and  515  (Figure  2.20)  maximized  and  for  the  blastward  and 
leeward  sidewalls  at  the  time  that  the  vertical  stresses  in  elements  369 
and  679  maximized  (i.e.,  when  bending  moment  maximized  at  the  midsections 
of  these  structural  elements).  The  shear  and  moment  diagrams  for  the 
roof  and  floor  are  plotted  in  Figures  3.11  and  3.12,  respectively,  while 
those  for  the  blastward  and  leeward  sidewalls  are  plotted  in  Figure  3.13. 
The  times  for  which  the  moment  and  shear  diagrams  were  plotted  for  each 
member  are  shorn  in  the  appropriate  figure.  The  signs  for  the  moments, 
shears,  and  thrusts  are  those  of  a right-handed  sign  convention. 

In  Figure  3.11,  the  maximum  moment  occurs  near  the  center  of  the 
roof  of  the  structure  fo-  all  calculations  and  is  positive.  The  maxi- 
mum moment  in  the  floor  of  the  structure  is  negative  and  occurs  slightly 
downstream  from  the  center  of  the  floor  of  the  structure  as  shown  in 
Figure  3.12.  The  maximum  moments  in  the  sidewalls  occur  near  the  top 
of  the  sidewalls  as  shown  in  Figure  3.13.  The  general  trend  is  for 
moments  to  increase  with  increase  in  weapon  yield  and  with  loose  backfill. 
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Hie  calculated  moments  and  thrustB  for  specific  sections  are  compared 
with  thrust-moment  interaction  diagrams  in  Figures  3. 11*  and  3.15.  Fig- 
ure 3.1*+  shows  a thrust-moment  interaction  diagram  for  yield  of  a struc- 
tural section  that  might  be  typical  of  those  for  the  roof  and  floor  in 
the  idealized  structure,  and  Figure  3.15  presents  the  same  type  thrust- 
moment  interaction  diagram  for  the  sidewalls.  The  thrust-moment  inter- 
action diagrams  were  determined  by  the  methods  of  Section  8.3  of  ASCE 

Q 

Manual  427  and  by  assuming  a concrete  strength  of  4.0  ksi  (277  bars),  a 
yield  strength  of  steel  of  52  ksi  (3600  bars),  and  1 percent  of  rein- 
forcing steel  in  each  face.  Also  shown  in  Figure  3.14  are  the  maximum 
thrusts  and  moments  on  Sections  FF  (through  the  roof)  and  Section  RR 
(through  the  floor),  and  in  Figure  3.15 » those  on  Sections  wu  and  KK 
(through  the  sidewalls).  As  further  noted  in  Figure  3.14,  an  increase 
in  weapon  yield  and  a change  from  dense  to  loose  backfill  caused  an  in- 
crease in  thrusts  and  moments.  The  same  trend  was  seen  in  Figure  3.15 
with  the  exception  of  Section  KK  for  the  10-KT  case.  The  increases  in 
thrust  and  moments  due  to  changes  in  the  type  of  backfill  were  larger 
for  1-MT  cases  than  for  the  10-KT  cases. 

What  do  the  increases  in  thrusts,  moments,  and  shears,  which  have 
been  shown  to  occur  as  backfill  quality  decreases,  mean  in  terms  of 
structural  integrity,  strength  design  criteria,  and  resulting  costs? 

The  structure  idealized  here  would  probably  undergo  significant  perma- 
nent deformation  and  might  even  be  totally  unusable  if  subjected  to 
the  loading  assumed.  Because  the  idealized  structure  is  linear  elastic, 
the  calculation  cannot  possibly  demonstrate  that  type  of  behavior.  Hie 
trends  in  Figure  3.14  indicate  that  structural  integrity  could  possibly 
be  threatened  by  an  increase  in  weapon  yield  or  decreases  in  backfill 
quality.  (This  trend  is  more  evident  for  the  roof  than  for  the  floor 
of  the  structure.)  The  some  statement  can  be  made  concerning  the  trends 
shown  in  Figure  3.15.  (For  the  sidewalls,  the  trend  is  more  evident  for 
the  1-MT  loading  ca3e.) 


3.2  SELECTED  TIME  HISTOKIES 

Figures  3.6  through  3. ^2  present  comparisons  of  selected  stress  or 
motion  time  histories  for  the  four  calculations.  In  each  figure,  the 
element  and/or  node  locations  are  sketched,  'lhe  sign  convention  for 
these  figures  is  as  follows:  upward  and  outward  movements  and  tensile 
stresses  are  considered  positive. 

Comparisons  of  the  vertical  displacement  and  velocity  time  histories 
for  the  center  of  the  roof  span  are  shown  in  Figures  3.16  and  3. IT*  re- 
spectively. As  can  be  seen  in  the  figures,  both  peak  downward  velocities 
and  peak  downward  displacements  Increase  for  loose  backfill  conditions  and 
higher  weapon  yields.  The  maximum  displacements  varied  more  with  the 
weapon  yield  than  with  the  backfill  properties,  while  for  the  peak  parti- 
cle velocities,  the  reverse  was  true.  After  the  first  relative  maximum 
displacement,  the  period  of  the  oscillation  in  the  motion  of  this  point 
ranged  from  60  to  05  ms.  ‘Hi is  corresponds  to  a frequency  range  of  15  to 
17  Hz.  Within  or  very  close  to  this  same  range  are  the  first  (15.1)  and 
second  (16.2)  frequencies  determined  in  the  in  vacuo  modal  analyses*  and 
with  the  second  (16.7),  third  (17.2),  and  fourth  (17.^)  frequencies  de- 
termined during  the  embedded  nodal  analyses**  using  the  upper  bound  mate- 
rial properties.  The  deflected  shapes  of  the  structure  at  the  times  that 
the  verticaL  deflection  at  thiB  point  maximize  and  at  the  100  msec  time, 
as  shown  in  Figures  3.6  and  3.6,  respectively,  can  be  approximated  as  a 
combination  of  the  mode  shapes  associated  with  these  frequencies. 

Figures  3.16  and  3.19  show  the  horizontal  stresses  at  points  three- 
fourths  of  the  way  from  the  neutral  axis  to  the  top  and  bottom  fibers  of 
the  roof,  respectively.  As  can  be  seen,  the  peak  stresses  above  and  be- 
low the  neutral  axis  are  of  opposite  sign  and  are  notably  higher  for  the 
loose  backfill  and  higher  weapon-yield  cases. 

Hie  vertical  and  horizontal  velocity  time  histories  for  the  upper 
blastward  corner  of  the  structure  are  shown  in  Figures  3.20  and  3.21, 


* See  Figure  A. 3 of  Appendix  A. 

**  See  Figure  A. 7 of  Appendix  A. 
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respectively.  Uie  peak  downward  velocities  for  the  loose  backfill  cal- 
culations sire  slightly  higher  than  those  for  the  dense  backfill  calcu- 
lations for  both  weapon  yields.  The  peak  outward  velocity  for  the  10-KT 
loose  backfill  calculations  is  only  slightly  greater  than  that  for  the 
10-KT  dense  backfill  calculation;  however,  the  peak  outward  velocity  for 
the  1-MT  loose  backfill  calculation  is  much  greater  (a  factor  of  2 higher) 
than  that  for  the  1-MT  dense  backfill  calculation. 

Vertical  and  horizontal  displacement  time  histories  for  a point 
(node  3bl)  near  the  midheight  of  the  blast ward  sidewall  of  the  structure 
are  shown  in  Figures  3.22  and  3.23,  respectively,  'flie  type  of  backfill 
has  little  effect  on  the  vertical  displacements,  as  shown  in  Figure  3.22; 
however,  the  peak  downward  displacements  for  the  calculations  with  1-MT 
loading  are  approximately  50  percent  higher  than  those  for  the  calcula- 
tions with  10-KT  loading.  The  peak  inward  displacement  at  the  midheight 
of  the  structure  is  highest  for  the  loose  backfill  and  largest  weapon 
yield. 

Figures  3.24  through  3.29  compare  vertical  deflection  time  histories 
of  points  on  the  blastward  face  of  the  structuz*e  with  those  for  adjacent 
points  in  the  backfill  for  three  different  depths.  As  also  noted,  the 
maximum  vertical  deflections  of  points  on  the  structure  are  not  affected 
by  changes  in  backfill  or  vertical  location  along  the  wall.  (The  de- 
flections of  the  structure  for  the  calculations  with  10-KT  loading  and 
1-MT  loading  are  5.0  and  7.5  centimetres,  respectively.)  The  reason  for 
this  is  that  the  vertical  deflection  of  the  structure  as  a whole  is 
mainly  controlled  by  the  applied  loading  and  the  properties  of  the  rela- 
tively incompressible  material  upon  which  it  rests.  The  figures  further 
show  that  the  relative  displacements  between  the  loose  backfill  and  the 
structure  are  much  greater  than  the  relative  displacements  between  the 
dense  backfill  and  the  structure.  The  points  in  the  backfill  examined 
were  only  JO  feet  (3.05  metres)  from  the  structure,  and  backfill  deflec- 
tions in  these  figures  at  these  locations  are  probably  less  than  they 
should  be  because  the  backfill  elements  are  welded  to  the  structural 
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elements  at  the  soil-structure  interface  (Bee  Section  3.1.1).  Neverthe- 
less, the  important  influence  of  backfill  properties  on  relative  dis- 
placement sensitive  items,  such  as  pipes,  cables,  and  connections  to  the 
structure,  is  clearly  evident. 

At  middepth  of  the  structure,  shown  in  Figures  3.2b  and  3.27,  the 
structure  moves  down  with  respect  to  the  backfill  at  early  times.  The 
phenomenon  occurs  because  a significant  stress  has  not  propagated  to  this 
depth  through  the  backfill  materials,  while  significant  stress  has  been 
applied  to  the  roof  of  the  structure,  which  is  only  5 feet  (1.52  metres) 
below  the  ground  surface  and  is  very  quickly  transmitted  into  the  founda- 
tion causing  it  to  deform.  Hie  fact  that  for  certain  cases  the  structure 
has  been  shown  to  deflect  downward  with  respect  to  the  backfill  at  early 
times  could  be  a very  significant  consideration  when  designing  cable  or 
pipe  connections  to  protective  structures.  The  loose  backfill  and  higher 
weapon-yield  cases  present  the  most  severe  conditions  for  relative  displace- 
ment of  the  structure-backfill  system  in  either  vertical  directions. 

In  Appendix  B,  the  vertical  deflection  of  the  base  of  the  structure 
(the  average  of  the  two  bottom  corners)  is  compared  to  the  rigid  body 
motion  of  the  structure  as  predicted  by  methods  outlined  in  Reference  5 
for  the  l-MT/300-psi  (20.7-bar)  dense  backfill  case.  Good  agreement  is 
demonstrated,  and  the  structure  is  shown  to  be  deflected  downward  nearly 
twice  the  free-field  deflection  at  a point  that  is  the  same  depth  as,  but 
upstream  of  the  base  of,  the  structure. 

Figures  3.30  and  3.31  show  the  time  histories  of  horizontal  stress 
on  the  blastward  sidewall  at  points  near  its  top  and  midpoint,  respect- 
ively. The  horizontal  stress  arrives  sooner  in  the  dense  backfill  cal- 
culations than  in  the  loose  backfill  calculations  due  to  the  higher  wave 
speed  in  the  dense  backfill.  Both  the  peak  stress  and  the  impulse  per 
unit  area  delivered  to  the  structure  are  higher  in  the  loose  backfill 
cases. 

Figures  3.32  through  3.37  present  comparisons  of  selected  stress 
time  histories  of  stress  components  within  the  blastward  sidewall  for 
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all  the  calculations.  Figure  3.38  compares  vertical  stress  time  histories 
for  an  element  on  the  inside  of  the  leeward  sidewall  near  the  midheight. 
Figure  3.39  compares  horizontal  stress  time  histories  for  an  element  on 
the  top  of  the  floor  near  the  center.  In  all  cases,  the  maximum  stresses 
increase  with  an  increase  in  weapon  yield  and  for  the  loose  backfill  con- 
dition. 

Examination  of  the  time  histories  of  stress  and  motion  shows  the 
presence  of  frequencies  up  to  250  Hz  in  the  wave  forms  for  points  in 
the  structure.  The  highest  frequencies  that  can  possibly  be  faithfully 
represented  are  given  by  the  reciprocal  of  the  rise  time,  i.e.,  1/10  ms 
or  100  Hz.  As  a practical  matter,  however,  credible  frequencies  will 
probably  not  exceed  one-half  this  value  or  50  Hz.  The  FE  grid  for  the 
backfill  material  immediately  surrounding  the  structure  can  transmit 
frequencies  through  the  dense  backfill  up  to  about  70  Hz.  However,  for 
the  case  of  the  loose  backx'ill,  the  grid  can  only  transmit  frequencies 
up  to  about  25  Hz.  Thus,  the  cutoff  frequencies  are  on  the  order  of 
50  Hz  for  the  two  dense  backfill  calculations  and  25  Hz  for  the  two 
loose  backfill  calculations. 

3.3  SHOCK  SPECTRA  ANALYSES 

The  effects  of  backfill  variations  on  the  response  of  possible  equip- 
ment mounting  points  on  the  roof,  sidewalls,  and  on  the  floor  can  be 
considered  by  examining  the  shock  spectra  for  the  motion  wave  forms  cal- 
culated for  these  locations.  Figures  3.**3  through  3.**b  present  the 
2 percent  damped  shock  spectra  calculated  for  the  velocity  wave  forms 
in  Figures  3.17  and  3.14  through  3.^2 . The  inside  of  the  roof  of  the 
structure  (node  513)  has  the  moBt  severe  calculated  shock  environment  of 
those  points  examined.  In  the  credible  frequency  range,  the  shock  environ- 
ment worsens  for  the  loose  backfill  and  the  higher  weapon  yield.  The 
center  of  the  floor  of  the  structure  has  the  smallest  shock  environment 
(Figure  3.^6).  This  is  probably  due  to  the  fact  that  the  floor  rests 
on  the  undisturbed  stiff  shale  and  is  somewhat  shielded  from  the  direct 
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effects  of  airblast-induced  ground  shock. 

The  maximum  response  of  the  sidewalls  and  the  roof  of  the  structure 
over  the  credible  frequency  range  occurs  at  frequencies  of  15  to  IT  Hz. 
This  frequency  range,  of  course,  coincides  with  the  frequency  range  ob- 
served for  structural  motion  time  histories  and  supports  the  discussion 
concerning  mode  shapes  and  frequencies  in  paragraph  2 of  Section  3.2. 

The  maximum  response  of  the  floor  of  the  structure  occurs  at  frequencies 
ranging  from  12  to  15  Hz.  This  response  could  possibly  be  controlled  by 
the  properties  of  the  clay  shale  upon  which  the  floor  rests. 

3.4  SUMMARY 

Tne  results  of  the  calculations  presented  herein  demonstrate  that 
both  loose  backfill  and  higher  weapon  yield  tend  tc  increase  the  load  on, 
the  stresses  and  bending  moments  within,  and  the  deflections  of  the  struc- 
ture. 

The  results  of  the  shock  spectra  analyses  also  show  that  the  shock 
environments  for  the  roof,  sidewalls,  and  floor  of  the  structure  are 
worsened  by  the  loose  backfill  and  the  higher  yield  conditions. 

lhe  frequencies  observed  for  selected  motion  time  histories  of  the 
structure  and  the  deflected  shape  of  the  structure  at  selected  times 
indicate  that  the  structure  is  responding  as  in  the  first  two  in  vacuo 
modes  and/or  the  second,  third,  and  fourth  embedded  modes  as  calculated 
with  upper  bound  properties.  These  observed  frequencies  range  from 
15  to  17  Hz  and  are  well  within  the  range  of  credible  frequencies  of  the 
calculations  as  controlled  by  the  airblast  rise  time  and  duration,  the 
FE  grid,  and  the  controlling  wave  speeds  of  the  materials. 
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Table  3.1  Maximum  and  Final  Rotations  of  the  Chord  Formed  by  the  Two 
Bottom  Corners  of  the  Structure;  300-psi  Nominal  Over- 
pressure, Dynamic  1-MT  and  10-KT  Calculations . 


Rotation 


Loose  Backfill 
rad  time,  ms 

1-MT  Calc ulations 


Dense  Backfill 
rad  time,  ms 


Counterclockwise  max 

0.00164 

55 

0.00131 

44 

Return  to  upright  position 

0 

70 

0 

100 

Clockwise  max 

0.00155 

145 

0.00107 

173 

Final  (clockwise) 

0.00098 

200 

0.00090 

200 

10-KT  Calculations 

Counterclockwise  max 

0.00123 

52 

0.0010T 

45 

Return  to  upright  position 

0 

65 

0 

82 

Clockwise  max 

0.00246 

200 

0.00246 

200 

Final  (clockwise) 

0.00246 

200 

0.00246 

200 

Table  3.2  Shear,  Thrust, and  Bending  Moment  for  Selected  Sections 
Through  the  Structure  at  a Time  of  62.5  ms  after  the 
start  of  the  1-MT  Dense  Backfill  Calculation . 


CTION 

SHF  A 9 FOKC E 

AXIAL  THRUST 

9t.N0 1 NG  H0MFNT 

N/M 

N/M 

N-M/M 

AA 

-• 102376E  *0  A 

-.419787t*07 

«.473*55F*07 

-.8422096*07 

-«  2o0  G66fc*07 

.6160  23E*07 

CC 

-.558145E+07 

ina7*9t>07 

• 141442E*08 

on 

-• 233  285E  *07 

-.582594E*0b 

• 1907836*88 

EE 

-.U70Q6E  + 06 

-.404017E*06 

.199071E+08 

FF 

.16b*31t*97 

-.4237254*06 

,1916b6F*08 

GG 

. 339078E  *0  7 

-«  4368u<*t*0b 

.17t9??E*08 

HH 

• 59O'J60E*Q7 

114  743  F 707 

• 1 25  77c£*  0 8 

II 

• 877  034F  + 0 7 

29841 4F*07 

.mF589?£*07 

JJ 

. 1136 01E*C8 

-. 4431996*07 

-.67  2^  264*0  7 

KK 

• 27C  986E*G7 

-.197*736*08 

-.If  52  22r*  0 B 

LL 

.134079E  +97 

20*990E*08 

-• 7t()289F*Q7 

MM 

• 881771*  F + 3o 

-.213870F708 

» • bt-0  369F*0  7 

NN 

. 884999F ♦ ’o 

-.213b'>2t*0<} 

-,4orw.  M8F  + 07 

00 

-•8327916*07 

.1953216*07 

-.2?93C8F*'*7 

PP 

-•29bui9F*37 

• 280952E  *07 

-. 798783t*P7 

QQ  - 

-.10S5Q4E+07 

• 3272316*07 

-. 993bQPE*07 

RR 

• 32(  70 1F*«  6 

.35452 9F*07 

-. If  48116*08 

SS 

.119  ,'94E*07 

. 3457894*07 

-.  1PJ973E*08 

TT 

. 140339F*07 

. 34390 7t *87 

-• 9259l5c*07 

UU 

.1996<*9E*07 

.34910 3£*07 

-.818*  o4l* C 7 

YV 

.2600J4t*C7 

.357457E*07 

-.51 65  28t*  07 

MM 

.3"4686F*"7 

.3421*76*07 

- . 1°3J  bbE*07 

XX 

.68F24lL*07 

.155h8UE*07 

. ?r9488E*07 

YY 

• bO 16b8F+36 

209155F *8  8 

.6  <3323E*C7 

Z7 

-.bb9890F*0b 

-•  196btl9F*0  8 

.*>431f7F*r7 

PH 

-.978  J72F*0fi 

190  817t*0  8 

. 6 298  C 3E*  07 

OM 

-• 19971 SF  *0  7 

17  2 1*1  26*0  8 

. 9464«*3F*07 

Table  3.3  Shear,  Thrust,  and  Bending  Moment  for  Selected  Sections 
Through  the  Structure  at  a Time  of  67.5  ms  after  the 
start  of  the  1-MT  Dense  Backfill  Calculation. 


SECTION 

SH£A K FORCE 

A*IAL  THRUST 

BLN 0 INS  MOMENT 

N/M 

N/M 

N-M/M 

AA 

-.1054Q8E+08 

-. 37432 Tt  + 07 

-.4458416+07 

B9 

-.83C  .87E+07 

-.2Z9356u+C7 

• P28365E*r7 

CC 

-.54lb04F+07 

-• 250  7l6t  + 07 

• 140169E+08 

00 

-• 238525E+07 

-•265  022E+07 

• 1 84372E+P  8 

EE 

-.927898E+06 

-.266036E+07 

• 192P  59E  + P 8 

FF 

•93142GE+06 

-.245843E+07 

• 1894  986+68 

G6 

• 2813094+97 

-• 233132E+07 

. l-6633t  + «8 

MM 

• 574132E  + 07 

-.245642E+07 

. I274?7t  + f 8 

II 

• 8U96  2L*0r 

-•  29b£?5E+07 

• 4970  1UE+P7 

JJ 

• 11U148E  + 08 

-• 4 16C2IE+07 

-•57476PE+07 

KK 

• 229079E+07 

-•18413  93 

108205F+08 

LL 

• 16  3122t  + 87 

-.192156^^08 

-• 773167t+07 

MM 

• 12k  7u6£+07 

-. 19551 4F  + 08 

-.641475c+07 

NN 

.48738 l£+08 

-• 196719E+08 

-.474711E+97 

00 

-.828929E+07 

• 151226F+07 

-.24047Sc.  + ?7 

pp 

-.284u95E+07 

•233U62E+07 

-•792474F407 

00 

-.939922E+L6 

• 196075E+O7 

-• 989753t+P7 

RR 

•40C944L+05 

• 218987E+Q7 

- • 1P15  98E  + 0 8 

ss 

•619629E+06 

• 21273  3t +07 

-•994206F+07 

TT 

• 832944F+0o 

• 236858t+07 

-.921564E+07 

UU 

• 15Fi<73£  + u7 

• 275304E  + 07 

- • 8243C4E+07 

VV 

• 21 7>8  7F+07 

• 37  330  3E  + 07 

- • 561 876F+0  7 

MM 

.27CB79*+07 

• 424453E+Q7 

-. 229156E+97 

XX 

•696349c+97 

. 26t  22  6E  + C7 

• 174157r+f t 

VY 

• 355575F+05 

- • 21bb97t  + 0 8 

• 53166  *F  + 07 

Z? 

• • 186319E  + 06 

-.  20910ot.+  0 8 

•566889C+07 

PM 

-.bHP24CE+06 

24177 8L  + 08 

•b316C5£+07 

OM 

-.182424E+J7 

-.1849294+08 

• «>3J553F  + P7 

.1 
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Table  3.U  Shear,  Thrust, and  Bending  Moment  for  Selected  Sections 
Through  the  Structure  at  a Time  of  65  ms  after  the 
start  of  the  1-KT  Loose  Backfill  Calculation . 


SECTION 

SHEAR  FORCE 

AXIAL  THRUST 

9EN0INS  NONENT 

N/M 

N/M 

N-M/M 

AA 

-.141963E*Q8 

-•§20  878E*§7 

-•5o§0»6E*07 

89 

-•  108962E*3  9 

-•414780E+I7 

• 830  392E  + 07 

cc 

-•728103E+Q7 

•.276F00L+07 

• 183627E+88 

00 

• • 343906F *07 

-• 176086F  + 07 

• 241612E+0  8 

EE 

••136 180E+07 

-• 159918E+07 

• 254624E+0  8 

FF 

•932809E+06 

16172 8E+07 

• 267662E+Q8 

66 

• 329538E+07 

-• 1§4365E»07 

• 2*‘5845F^08 

HH 

.617370E+07 

-• 106 12  4E+07 

• 1 8**952E*Q  6 

II 

• 895037F  *07 

-•<*fl4154E*06 

• 9f  8375F+07 

J 1 

.128761E+38 

-• 979693t>05 

-• 383827E+07 

KK 

,560  826E+06 

-• 184591E+Q8 

1K529GE*08 

LL 

. 293360E+37 

-.177870E+08 

-•  130873E+08 

NM 

,339794E»07 

-.178012E*08 

107252E*08 

NN 

• 441 126E*17 

-• 186154F*08 

-.4flPC6E*07 

00 

-• 677665F+07 

• 747343E«-07 

•• 5397  21E+07 

PO 

•• 236499E*07 

•709?fllEt07 

-•  1C0923E+Q8 

OQ 

• • 65566 1E+06 

• 60952  3E*87 

-•121351t*Q8 

RP 

•453260E+06 

•61442SE+07 

-• 1210  88E+08 

SS 

• 1 12561E+07 

•589591E*07 

-.119727E408 

TT 

• 220010E+07 

•618898E+07 

-•1P9972E+08 

UU 

• 294998E+07 

•634900E+07 

-• 958381E* 07 

*v 

• 31Q014E+Q7 

•678928E+07 

-.581641E+07 

HH 

. 337758E+07 

• 627631E+07 

•• 2188 16E*07 

XX 

.444350E+07 

• 522094fc>07 

. 11450ft.  *07 

VY 

-.244188E+07 

-• 189039E*08 

•4f 2251E*07 

Z7 

-• 18 180 1F*0  7 

-• 200  37 1E*B8 

•7b2922E*07 

PH 

•• 229627E+37 

-• 205974E*I8 

• 893792E+07 

OH 

-.364»321E*07 

-• 20  4530E+0  8 

• 13«5893E*0  8 

* 
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Table  3.5  Shear,  Thrust,  and  Bending  Moment  for  Selected  Sections 
Through  the  Structure  at  a Time  of  JO  ms  after  the 
start  of  the  1-MT  Loose  Backfill  Calculation  , 


SECTION 

SHEAR  FORCE 

AXIAL  THRUST 

ENDING  MOMENT 

N/M 

N/M 

N-M/M 

AA 

-.1261006*08' 

94269 4E487 

-.3618966407 

89 

-.99  8 3406  * 07 

".463273E+07 

.9003636*07 

CC 

“•669  0106*07 

“• 331 368E+07 

• 1958  88E*0  8 

00 

-.286998E*07 

-.2638616*07 

• 2541 94E*08 

EE 

-.8547666*06 

-•244692E407 

.2662796*08 

FF 

• 1066636*07 

-•259695F407 

• 2688 19E  + 08 

66 

•3470186*07 

-,306256E*Q7 

.258397E*08 

HH 

•882408E*07 

-.4128746*07 

• 1918816*08 

II 

• 1209436*08 

-•*5167946407 

.6926816*07 

JJ 

.1496946*08 

-.6022646407 

-.8815516*07 

KK 

.4698366*07 

-•227296E408 

-.152649E*08 

LL 

.2866516*07 

-.236  2556*0  8 

-.9287016*07 

HM 

.108774F*07 

-.2353796408 

-.7750576*07 

NN 

.6995066*06 

-.2420166408 

- . 7C01 14E  + 07 

00 

-• 1 09216E  40  A 

• 26131 4E407 

- * 2958 17E*  07 

PP 

— .400  30 46*37 

.3706636407 

-•1C  23176*6  8 

00 

-.1131496*07 

•446511F407 

-. 1238036*08 

R9 

.36  5 9 516  * 06 

•401 £896407 

-. 125332E+08 

SS 

• 9 16 139E  406 

•3986546407 

-.1193396*08 

TT 

• 673871F.*06 

.4136706*07 

-.1127686*08 

uu 

•120Q50F487 

.4409776407 

-. 1C5219E+08 

VV 

• 2013416487 

.5413926407 

-.7778646*07 

ww 

•234823F+07 

•569884E407 

- • 484150E+07 

XX 

•561569E407 

•5741916*07 

-.3948626*06 

vy 

-.327899E+07 

-.2167496*08 

.534396F*07 

77 

-.2076016407 

-.2048126*08 

. 9704506*07 

PH 

-•223'»6A£407 

-.1970056*08 

. 11524bE*0  8 

OH 

-•1374326*0* 

-.1870486*08 

. 1442456*08 
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Table  3*6  Shear,  Thrust, and  Bending  Moment  for  Selected  Sections 
Through  the  Structure  at  a Time  of  55  ms  after  the 


start  of  the  10-JCT 

Dense  Backfill 

Calculation . 

SECTION 

SHEAR  FORCE 

AXIAL  THRLST 

BENDING  8CMEN7 

N/M 

N/M 

N-M/M 

AA 

-.7P8742E+07 

31«741E*<  7 

•♦29386'  €*C~ 

BB 

-.56C417E* 37 

*•  226F  E6E  *\)7 

♦ 398263E  *07 

CC 

“♦383C97E* J? 

19082-36*07 

. 9136516*07 

00 

-♦143C84E*07 

-.1732566*07 

• 1 1 92 13a*C  8 

EE 

-.344U9E*!>6 

-.  170868E*P7 

♦123111E*C8 

FF 

• 65  896  ?E  + 06 

196o43E*37 

♦ 12?437E*r  8 

GG 

.1  368596*17 

2147 58 E *0 7 

♦ 1158366*08 

HH 

♦ 25  98 1>  BE  *17 

-♦  2393546*87 

.9  372  12E*C>  7 

II 

♦524131E*37 

“,252498E*v7 

♦ S'  3r87E*»  7 

JJ 

♦823865E* 17 

-♦  31772  £*.:7 

-♦3.117'“E*C7 

KK 

♦133345E+17 

12f<  698  6*  ' 8 

“.bOJ'O'-FU? 

LL 

.127fc4fcE*J? 

-.1169146*38 

-.6  389  3?E*P7 

MN 

♦1 376Q4E*37 

1 13797  E 8 

“•5  434fc7E*07 

NN 

♦ 1 4 168  7F  + 17 

-.  1L.92276**'  a 

-.3  HIJElfr 

00 

- .4  3471 CF * J 7 

. 1635976*07 

• . 2 43  ^ 90  E*(i7 

PP 

-♦2P  88S5E*37 

• 18051 5£  *0  7 

-.6 1C922E+C? 

00 

-.iri039E*J7 

♦ 2441676  *0  7 

-.7  891 29E*t  7 

RR 

.41T891E06 

. 295934EMJ7 

-.8  ?521'?E*f  7 

ss 

♦ IP  177 4E +97 

♦ 295137E  *07 

-♦8  « C 27?E*( 7 

TT 

♦169P46E+17 

. 312223  £+0  7 

-,7  i955qf > 

UU 

♦2C4098Z+J7 

. 31 742 5 E +7-7 

- ,i  8 88  24£  ♦!  7 

VV 

.2691396*17 

♦3v0242E*17 

-.3125  72  E*f 7 

MM 

♦2742lCE*37 

♦ 256782E  *P  7 

-.985  27r  E*Ct 

XX 

• 3 86f  4 8E*07 

♦ 1C  77C IE  *0  7 

♦2611 176*07 

YY 

.30  298t>f +0 5 

“•1  26995c  ♦<}  * 

♦ 3 710  766*0  7 

ZZ 

3713556*06 

“.125t8Tfc*-*  8 

♦ 4 f 38i«6F*i  7 

PH 

“ ♦ 15  715  2E*  J6 

12428 ‘ t *t-  8 

♦435202F*C7 

OH 

“.114t92E*37 

“ ♦ 1197516*08 

.5247576*17 

Table  3.7  Shear,  Thrust, and  Bending  Moment  for  Selected  Sections 
Through  the  Structure  at  a Time  of  60  ms  after  the 
start  of  the  10-KT  Dense  Backfill  Calculation . 


SECTION 

SHEAR  FORCE 

axial  thrust 

3ENDING  MOMENT 

N/M 

N/M 

N-M/M 

AA 

-.667139E*'7 

-.2?F$G2E4$7 

-.279653E+C-7 

89 

-.5?9*i2E+07 

-.132097E+07 

.445817E+1  7 

CC 

-.368965E407 

-*f 81998E+06 

» 1 OG  22'  F + C 8 

00 

-.129789E+J7 

.1858S1C4"* 

.1  29€8tE4t.6 

EE 

•22E71 8E+ lb 

• 1" 2783E4C6 

.13384'  E4C3 

FF 

• 1 47794E4Q7 

• 482639E  *0  v 

• 1 2983?  E 4G  8 

GG 

.268*16E417 

*44f 42?E**5 

• 1 178  32E  +06 

HH 

.Lbfir  36E*J7 

• • 295944E  4f  6 

.752211E4T  7 

II 

•59821 9E *17 

-.155t*.8F4f> 7 

.177484E4C  7 

JJ 

.52479LE4I7 

-•2S9b86E  *0  7 

-.493959L4C7 

KK 

• 129277E07 

-•11782?E*0S 

-.654548E*C7 

LI 

• 461C67E+J6 

• • 1 23938F  8 

-.511871E4T 7 

HN 

.1  l^rsF^ie 

-•1 2713 1 E+.^S 

-*474€34E4f 7 

NN 

*35l49EE4)fc 

••129 67 4 £4.'  f» 

• • 4 P52  6EF  4L  7 

00 

-.56i722E+37 

• b 78942  E *0  t 

-*17641uF4C7 

PP 

- • 2 44739c* J7 

. 1 1 486bc  +U 7 

-.621389E4T7 

QQ 

-.952297E+16 

. 152ll5E4r  7 

• • 7 9844F  E+f  7 

RR 

.138497E+D6 

• 19 i 272E+07 

-.842318L4L7 

SS 

•52E538E+16 

. ?t  9664E  4fl  7 

-.36167FE447 

TT 

.1C  V 66E  + 37 

• 231309E  4P7 

-.743547F*<  7 

UU 

.If f 731;*37 

• 2C  v*7ti4c  40  7 

-.639618E41  7 

VV 

•21C882E*37 

• 3 » 3057E  4f  7 

-.35987*  E*'7 

MM 

• 2 H98fl 1E  + J 7 

. 23863SE407 

”»E7320QF4Ct 

XX 

.431748E+37 

• 924269  E *0  t 

.274iaeE4C7 

YY 

.729793E+1o 

• • 1 42492E  *G  8 

*3  965  7GF4f  7 

ZZ 

-.247f 14E+ifi 

14363?E4**e 

.183;'  5r*E4<  7 

PH 

• *5  C 8 11  ?E *‘7 6 

-.1 1584tF*C8 

.436164E4C7 

OM 

-.iniritE+37 

U8931£*C<- 

• 5 893o!E  + C ~ 

Table  3.8  Shear,  Thrust, and  Bending  Moment  for  Selected  Sections 
Through  the  Structure  at  a Time  of  60  ms  after  the 
start  of  the  10-KT  Loose  Backfill  Calculation  . 


SECTION 

SHEAR  FORCE 

AXIAL  THRUST 

RENOING  MOMENT 

N/M 

N/M 

N-M/M 

AA 

-.8370636+07 

-.251423E+07 

-•401 623E+  07 

8* 

-.681801E+07 

-• 1523481 +07 

•4775276+07 

CC 

-.5059236+07 

-.5480726+06 

• 115312E+08 

00 

-.267273E+07 

• 3905366+06 

. 1588  466  + 08 

E£ 

-.1288646  +07 

.6581176+06 

.1693976+08 

FF 

. 876684E+05 

• 644700  F.  + 0 6 

• 1727  88E  + 0 9 

GG 

• 1477246+07 

.6009366+06 

• 167531E+0" 

HH 

.4564436+07 

• 38220  56  + 06 

• 1295576+0  8 

II 

.6591446+07 

.6166136+05 

.5769936+07 

JJ 

•7P2297E+37 

-.4722936+06 

27o2*4E+C7 

KK 

••761 Q95t +06 

-.QC2724L+07 

-.9*2G  87E  + C7 

LL 

• 76587 8E+0b 

10020  9E+08 

-•1009276+08 

MM 

• 260  108E+07 

-. 105343E+08 

-. 885179E+07 

NN 

•4179466 ♦07 

-• 115F81E  + 0 8 

-.  3(04206  + 07 

00 

387?88F+07 

• 59177  IF  + 07 

-.  4-'J627F  + Q7 

po 

-.1577526+97 

.6107696+07 

-• 8 C8221E  + C7 

QQ 

•200900E+06 

•6245T0E+07 

-.9246546+07 

RP 

• 1 31 8 10E+0  7 

. 61410 2E+07 

-.  8..8592E  + 07 

SS 

• 154511E+07 

•583787E+07 

-• 779471E+07 

TT 

• 19123OE+07 

. 54464 5E+07 

-.6Q2124E+07 

uu 

• 190 1256+07 

•6014536+07 

-.5860366  + 07 

YV 

.211557E+07 

.4304226+07 

-.3395836+07 

WM 

• 292207E+07 

• 294156E+07 

-. 1P3812E+07 

Xx 

.37C098E+Q7 

• 18501 6E+97 

. 2121 83E+U7 

VY 

- • 4756726  *06 

-.1249286+08 

• 3304156+07 

ZZ 

-.125977E+C7 

-.12193JE+08 

.4524736+47 

PH 

••  1577Q6E+07 

-•1218916+08 

.5574286+07 

OM 

-•194803F+07 

-.1159756+08 

. 845751E+07 
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Table  3.9  Shear,  Thrust, and  Bending  Moment  for  Selected  Sections 
Through  the  Structure  at  a Time  of  62.5  nis  after  the 
start  of  the  10- KT  Loose  Backfill  Calculation  . 


It* 

■ 


SECTION 

SHFAR  FORCE 
N/M 

AXIAL  THRUST 
N/M 

9ENDING  MOMENT 

N-M/M 

AA 

-.930198E+07 

-• 347 17  3t+Q7 

-.3672526*07 

89 

-.719139E+07 

-.300 1886*07 

• 603143E+07 

cc 

-.477745E+07 

-• 230685E*07 

. 1279016*08 

00 

-.2Q7519E+07 

-.1620046*07 

. 1656626*08 

EE 

-.728168E+06 

-.1253746*07 

. 1725026*0  8 

FP 

• 79931 3E  + 96 

-• 9850636*0  6 

• l7137KE  + n 8 

GG 

•192589E+07 

-.64530 16*06 

.1613676*08 

HH 

.39C547P+07 

-.300  3096*06 

• l?4u  87E+0  8 

IT 

•61Q161E+07 

-. 1807326*06 

. 6332136*07 

J J 

.884546E+07 

225686E+06 

3 **53966*07 

KK 

•5421276*06 

-• 1198226*08 

-.1045746+08 

LL 

• 154763E  *07 

-. 116976E+B8 

-.8840296*07 

MM 

.2261666*07 

-• 113745E+08 

« . 755339c*0  7 

NN 

• 299993E+97 

-.lu8358L*0fl 

-.36+0  68E+07 

00 

-.4287386*07 

•3807116*07 

- . 3510  84E+D7 

PP 

-.1712376*07 

.40  36376*07 

-.73o229E*07 

00 

-.84565 76*08 

. 42493 8E*07 

-.  894381E*07 

RR 

* 250  787E+05 

•4863146*07 

- • 9380  42E  + U 7 

SS 

• 923567E*06 

•515836E*07 

-.893055^*07 

TT 

.1786876*07 

• 5184656*07 

-. 8204256*07 

uu 

. 2452626*07 

,530  2286*07 

- . 666940E+07 

.23<U31E*07 

• 474244E*0  7 

-• 4112  92E  + G 7 

VM 

.299119E+P7 

• 3736106*07 

-.  1656  3BE+0 Y 

xy 

• 302t  89E*07 

• 3319896*07 

. 2113666*07 

YY 

-•  147  808E*07 

-• 1Z6736E*0  8 

. 3474806*07 

Z 7 

-•  106456F+P  7 

-•  13u 2856*0  8 

• 54j 1366*07 

PH 

-.1455096*07 

-• 133681E+08 

. 6262156*07 

OH 

-•1947116*07 

-• 129644E+08 

. 90Q207E*07  * 

68 


Table  3*10  Shear,  Thrust, and  Bending  Moment  for  Selected  Sections 
Through  the  Structure  at  a Time  of  67 «5  ms  after  the 
start  of  the  10-KT  Loose  Backfill  Calculation  . 


SECTION 

SHEAR  FORCE 

AXIAL  THRUST 

9EN0ING  MOMENT 

N/M 

N/M 

N-M/M 

AA 

-.774121F  »07 

-.292314E+0? 

-416717*E<*07 

B1 

-.5<*6176F*07 

-• 293576E  *07 

•6f6854E*87 

CC 

-• 382717E*07 

-. 2614535*07 

. 1258725*08 

00 

-.137209E*07 

-• 217  A55£*07 

• 1&6E  33E*0  * 

EE 

- • 1U8191E+6 6 

-.2271785*87 

. 1635425*98 

FF 

•131357E407 

2511085*87 

• 16?448E*06 

GG 

.281744E*07 

-.2712lOE*87 

•1?1872E*08 

HH 

.512011E*07 

-.321084E*07 

. 1"9112E*0* 

II 

«724939F*07 

-. 37735 7E*87 

• 3750 25E*07 

J J 

.838037E  + 07 

-.40  2683E  *07 

• * 5*>0  3795*07 

KK 

.241597E407 

-• 1274745*08 

-.882071E*07 

LL 

•120910E407 

-• 1367235*08 

5°5524E*07 

NM 

•404682E+06 

14265 JF*08 

-•53130  7E*tl7 

NN 

• 19A027E4|)6 

-.1&0613E+OA 

-.515472E*07 

00 

- • 672730E  *0  7 

• 3u4 10  It 406 

-.194052E*07 

PP 

-.27fl490E*07 

• 563297E  +0f> 

-. 690150fc*P7 

00 

-.895387E406 

. 947876E+06 

-.873054E*07 

RR 

-.728234E+05 

•113034t*07 

-.9C73?5E*07 

ss 

.44822?E*06 

. 12141 5F*87 

-.884865E*P7 

TT 

• 910445EO6 

. 14350 1£*07 

-.  83*.929E*07 

uu 

. 134963E+07 

• 17721 15*07 

-.759R09E*07 

vv 

.1!>7545E»07 

• 322431E  407 

-.5H75C5  + P7 

WM 

• 291679Ft07 

.337939E*07 

-.3160765*07 

xy 

• 31249UE  + 0 7 

.4275645*07 

. 1205675*07 

vv 

- • 2P8479E  *0  7 

-• 1 322605*03 

• 3f  694  AE*07 

71 

-.134446F+07 

129743E*B8 

•6b4*52E*07 

PH 

-.906459E*06 

-.127?76E*08 

• 7b4057t*07 

OH 

-.118 194F4A6 

-• 119210F*08 

• 837*1  61F*07 

4 


4 


j 


: 

I 


69 


4 


Jm 


■ 


gur< 


NOTE  I MT  LOOSC  no  DEFLECT-ON  AT  T.ME 


DC  FLECTION 


scale 


o.  DEFLECTED  SHAPE  OF  STRUCTURE 


normal 
STRESS  , 


NORMAL 
STRESS  , 


normal 
STRESS  . 


b.  STRESSES  NORMAL  TO  STRUCTURE 

Figure  3.5  Comparison  of  deflected  shape  of  and  stresses 
normal  to  the  structure  at  Uo  ms  after  the  be- 
ginning of  the  problem. 
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b.  STRESSES  NORMAL  TO  STRUCTURE 

Figure  3*10  Comparison  of  the  deflected  shape  of  and 
stresses  normal  to  structure  at  200  ms 
after  the  beginning  of  the  problem. 
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n<»re  Comparison  of  maximum  moment  and  shear  diagrams  for 

selected  sidewall  sections  at  the  time  that  the 
vertical  stress  in  elements  369  and  6j 9 maximise  in 
the  blastvard  and  leeward  sidewalls,  respectively. 
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3.15  Moment-thrust  interaction  diagram  for  a rectangular  beam-column  and  data  for  selected 
sidewall  sections. 


YIELD  BACKFILL 


Fi^jure  3.16  Vertical  displacement  tine  histories  for  the  center  of  the  structure  roof. 


time  y msec 

Vertical  velocity  tine  histories  for  the  center  of  the  structure  roof 


LEGEND 


TIME,  MSEC 

Horizontal  stress  time  history  for  a point  in  the  roof  three-quarters  of  the 
way  from  the  neutral  axis  to  the  bottom  fiber  (element  5lU). 
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Vertical  velocity  tioe  histories  for  a point  (node  268)  on  the  blastvard  sidewall 


outward  LEGEND 


GROUND  SURFACE 


TIME,  MSEC 
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Fig Lire  3.39  Horizontal  stress  tine  histories  for  elerent  515  at  the  center  of  the  floor 
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FiS-^re  3.44  Vertical  shock  spectra  for  a point  (node  >13)  on  the  inside  of  the  structure 
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? igure  3.^5  Horizontal  shock,  spectra  for  a point  (node  666)  on  the  inside  of  the  leeward 
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CHAPTER  4 

PRESENTATION  AND  ANALYSES  OP  RESULTS  FOR  10-KT  WEAPON  YIELD, 

1000- PSI  OVERPRESSURE  DENSE  BACKFILL  CALCULATION 

4.1  CALCULATION  DESCRIPTION 

An  additional  dynamic  calculation  was  conducted  identical  to  those 
described  in  Chapter  2,  except  the  surface  boundary  of  the  problem  was 
loaded  with  a traveling  airblast  from  a 10-KT  weapon,  which  produced  a 
nominal  1000-psi  (68.9-bar)  overpressure  loading  on  the  center  of  the 
structure.  This  calculation  was  conducted  for  the  dense  backfill  case 
only.  The  applied  surface  overpressure  time  histories  over  the  upstream 
and  downstream  boundaries  of  the  grid  and  over  the  center  of  the  struc- 
ture are  shown  in  Figure  4.1.  There  is  a very  large  gradient  in  peak 
surface  overpressure  (lOx)  over  the  240-foot  (73.2-metre)  range  covered 
by  the  calculation  grid  as  indicated  in  the  figure.  No  artificial  rise 
time  was  appended  to  the  front  of  this  loading  function  because  the 
added  impulse  for  a 10-KT  yield,  1000-psi  (68.9-bar)  overpressure  rela- 
tion would  be  too  great  if  a rise  time  commensurate  with  the  frequency 
transmission  capability  of  the  grid  was  appended.  A finite  rise  time 
will  develop  in  the  media  due  to  the  grid  limitations  after  the  airblast 
has  propagated  into  it  for  some  small  distance.  Credible  frequencies 
probably  will  not  exceed  50  Hs.  The  maximum  objective  of  adding  this 
calculation  was  to  gain  qualitative  insight  into  the  effect  of  this 
additional  case  on  the  external  loading  on  the  structure. 

4.2  STRESS  AND  DISPLACEMENT  AT  SELECTED  TINES 

4.2.1  Deflections  Across  Backfill  Sections 

Calculated  vertical  displacement  patterns  for  sections  through  the 
dense  backfills  at  40  and  100  ms  are  shown  in  Figure  4.2.  At  the  5-foot 
(1.52-metre)  depth,  the  deflection  pattern  on  the  upstream  side  of  the 
structure  is  the  same  for  both  the  40-  and  the  100-ms  time.  However,  on 
the  downstream  side  of  the  structure,  the  backfill  deflections  at  100  ms 
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are  greater  than  at  40  ms  because  at  UO  msec  the  peak  stress  has  not  yet 
propagated  to  the  5-foot  (1.52-metre)  depth.  Similar  behavior  is  also 
shovn  for  the  15-foot  (4.57-metre)  depth.  The  maximum  backfill  deflec- 
tion at  the  5-foot  ( 1.52-metre)  and  15 -foot  (4. 5 7-metre)  depths  in  the 
vicinity  of  the  structure  at  100  ms  are  23  and  12  centimetres,  respect- 
ively. For  the  10-KT  weapon  yield,  nominal  300-psi  (20.7-bar)  loading 
over  dense  backfill  calculation,  the  maximum  deflections  at  the  same 
depths  and  time  were  19  and  7 centimetres,  respectively. 

4.2.2  Stresses  on  and  Deflection  of  the 
Exterior  Surfaces  of  the  Structure 

Instantaneous  distribution  of  normal  stresses  on  the  exterior  sur- 
faces of  the  structure  at  40,  70,  150,  and  200  ms  are  shown  in  Figure  4.3. 
Instantaneous  patterns  of  deflection  along  the  exterior  surfaces  of  the 
structure  at  the  same  times  are  shovn  in  the  same  figure.  The  40-ms 
time  is  close  to  the  time  (45  ms)  that  the  maximum  deflection  of  the 
center  of  the  roof  occurs.  At  40  ms,  the  stresses  are  concentrated  under 
the  sidewalls.  The  floor  and  the  roof  are  bowing  inward,  and  the  walls 
are  bowing  outward.  At  70  ms,  the  high  stresses  under  the  sidewalls 
have  been  relieved,  the  roof  is  rebounding,  and  the  downward  deflection 
of  the  center  of  the  floor  is  catching  up  with  that  of  the  sidewalls. 

The  blastvard  sidewall  is  undergoing  a shear  deformation,  while  the 
leeward  sidewall  has  about  rebounded  approximately  to  its  undeformed 
shape.  At  150  and  200  ms,  the  stress  on  the  structure  is  relieved,  and 
the  complete  structure  is  moving  upward,  taking  a shape  similar  to  its 
unstressed  shape  and  rotating  in  a clockwise  direction.  The  patterns 
of  calculated  stresses  on  and  the  deflected  shapes  of  the  exterior  sur- 
faces of  the  structure  observed  here  are  qualitatively  similar  to  those 
observed  for  the  structure  for  the  10-KT,  300-psi  (20.7-bar)  dense  back- 
fill calculation  whose  results  are  shown  in  Figures  3.6  through  3.10. 

4.2.3  Rotation  of  the  Structure 

The  structure  follows  the  same  trends  in  rotation  during  this  cal- 
culation, as  shovn  in  Table  4.1,  as  it  did  during  the  nominal  300-psi 
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(20.7-bar)  overpressure,  dense  backfill  calculation  (Table  3.1).  The 

structure  rotates  initially  in  a counterclockwise  direction  and  then  in 

a clockwise  direction.  The  maximum  rotation  in  this  calculation  was 

0.00328  radians  and  is  approximately  one-third  greater  than  the  maximum 

rotation  in  the  10-KT/300-psi  (20.7-bar ) /dense  backfill  calculation. 

4.2.4  Thrusts,  Shears,  and  Bending  Moments 
Within  the  Structure 

Axial  thrust,  shear,  and  bending  moment  per  unit  depth  of  the  struc- 
ture were  calculated  from  the  output  data  and  tabulated  at  10-ms  intervals 
for  the  28  structural  sections  shown  in  Figure  2.20.  From  the  data  in 
Table  4.2,  shear  and  moment  diagrams  were  plotted  for  the  roof,  floor, 
and  sidewalls  for  times  very  near  those  when  the  axial  stress  in  ele- 
ments near  the  center  of  each  member  maximised.  This  time  was  found  to 
be  approximately  40  ms  for  all  components.  Figure  4.4  presents  moment 
and  shear  diagrams  for  the  structure  floor  and  roof  at  40  ms,  while 
Figure  4.5  shows  those  for  the  sidewalls. 

Bie  trends  in  the  moment  and  shear  diagrams  for  the  roof  and  floor 
of  the  structure  are  the  same  for  this  calculation  (Figure  4.4)  as  they 
were  for  the  calculations  with  the  300-psi  (20. 7-bar)  nominal  overpressure 
loading  (Figures  3.11  and  3.12).  The  maximum  moment  at  the  center  of  the 
roof  from  the  300-psi  (20.7-bar),  10-KT  dense  backfill  calculation 
(Figure  3.11)  was  13  x 10b  Newton-metres/metres,  The  maxi  ..urn  moment  at 
the  center  of  the  roof  from  this  calculation  is  approximately  20  x 10b 
Newton-metres/metres.  Therefore,  the  effect  of  raising  the  applied 
overpressure  from  300  psi  (20.7  bars)  to  1000  psi  (68.9  bars)  for  this 
case  increased  the  moment  at  the  center  of  the  roof  by  a factor  of 
approximately  1.5.  However,  the  maximum  moment  for  this  case  was  not 
as  high  as  that  for  the  1-MT/ 300-psi  (20.7  bar) /loose  backfill  case  and 
was  about  equal  to  that  of  the  1-MT/ 300-psi  ( 20.7-bar ) /dense  backfill 
case.  The  maximum  moment  near  the  center  of  the  roof  occurred  at 
Section  EE  for  this  calculation,  and  at  Section  FF  for  the  previous  cal- 
culations. 
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The  trends  of  the  moment  and  shear  diagrams  for  the  structure  side- 
walls,  as  shown  in  Figure  4.5,  for  the  10-KT/1000-psi  (68.9-bar )/dense 
backfill  calculation,  are  similar  to  those  for  the  sidewalls  from  the 
previous  calculations  (Figure  3.13).  The  maximum  moments  in  the  blast- 
ward  and  leeward  sidewalls  for  this  calculation  occur  at  Sections  uko 
and  KK,  respectively,  were  9.5  x 10°  and  10.8  x 10°  Newton-metres/metres, 
respectively,  and  are  factors  of  1.6  and  1.3,  respectively,  greater  than 
those  from  the  10-KT/300-psi  ( 20. 7-bar) /dense  backfill  calculation.  Hie 
maximum  shear  in  the  leeward  sidewall  from  this  calculation  is  almost 
identical  to  that  from  the  l-MT/300-psi  (20.7-bar ) /dense  backfill  calcu- 
lation. 

The  maximum  moments  and  thrusts  that  occurred  in  the  structure  roof, 
floor,  and  sidewalls,  as  shown  in  Figures  4.4  and  4.5  and  in  Table  4.2, 
are  tabulated  below: 

Thrust 

N/m  (kip/ in.) 


Structural 

Section 

Moment 

Member 

(Figure  2.20) 

N-m/m  (kip- 

■in. /in 

Roof 

EE 

r 

20.7  x 10g 

(4650) 

Floor 

RR 

11.4  x 10? 

(2580) 

Blastward 

mm 

9.4  x 10° 

(2130) 

Sidewall 

f. 

Leeward 

KK 

10.8  x 10° 

(2420) 

0.17  x 10°  (1) 

5.2  x 10  ,(30) 
19.88  x 10  (115) 

17.9  x 10°  (100) 


Sidewall 


These  values  of  thrust  and  moment  are  very  close  to  those  for  the  1-MT/ 
300-psi  (20.7-bar ) /dense  backfill  calculation. 

4.3  SELECTED  TIME  HISTORIES 

Figures  4.6  through  4.11  present  comparisons  of  selected  stress  and 
motion  time  histories  from  the  calculation.  In  each  figure,  the  element 
and/or  node  locations  are  sketched.  The  sign  convention  for  these  fig- 
ures is  as  follows:  upward  and  outward  movements  and  tensile  stresses 
are  considered  positive. 

Figures  4.6  and  4.7  show  vertical  displacement  and  vertical  velocity 
time  histories,  respectively,  for  points  on  the  center  of  the  structure 
roof  and  floor.  The  maximum  downward  displacement  of  the  center  of  the 
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roof  (which  includes  some  rigid  body  displacement  of  the  structure)  is 
10. 4 centimetres  and  occurs  at  a time  of  45  ms  after  the  start  of  the 
problem.  The  peak,  downward  displacement  for  the  center  of  the  roof  from 
the  lO-KT/100-psi  (68.9  bars)/dense  backfill  calculation  (Figure  4.6)  is 
about  the  same  as  that  from  the  l-MT/300-psi  (20. 7-bar) /dense  backfill 
calculation  (Figure  3.1b).  The  peak  downward  velocity  for  the  center  of 
the  roof  is  0.77  centimetres/ms . This  value  is  approximately  50  percent 
higher  than  that  from  the  l-MT/300-psi  (20. 7-bar) /dense  backfill  calcu- 
lation. After  the  first  relative  maximum  in  displacement  occurs,  the 
period  of  oscillation  at  the  center  of  the  roof  (node  513)  was  approxi- 
mately 55  ms  (i.e.,  18  Hz).  Hiis  is  slightly  higher  than  the  frequency 
range  observed  for  the  previous  four  calculations  (see  Section  3.2  of 
Chapter  3). 

Figures  4.8  and  4.9  show  the  horizontal  displacement  and  velocity 
time  histories,  respectively,  for  points  at  the  center  of  the  sidewalls. 

The  trend  in  the  horizontal  displacement  versus  time  for  the  center  of 
the  blastward  sidewall  is  similar  to  that  shown  in  Figure  3.23  for  the 
l-MT/300-psi  ( 20 . 7-bar ) /dense  backfill  calculation.  The  center  of  the 
leeward  sidewall  is  always  deflected  leewardly  from  its  original  position. 

Figure  4.10  shows  the  horizontal  stress  time  histories  for  points 
near  the  inside  extreme  fibers  at  the  approximate  centers  of  the  structure 
roof  and  floor.  Both  the  points  are  in  tension  at  40  :na.  This  is  expected 
from  the  deflected  shape  of  the  structure  (Figure  4.3).  'Hie  maximum  ten- 
sile horizontal  stress  for  element  514,  near  the  center  of  the  roof,  is 
430  bars,  which  is  about  the  same  as  the  maximum  tensile  stress  at  the 
same  point  from  the  l-MT/300-psi  (20. 7-bar) /dense  backfill  calculation 
(Figure  3.19). 

Figure  4.11  shows  the  vertical  stress  time  histories  for  points  near 
the  inside  fiber  and  midheights  of  the  sidewalls.  The  maximum  compres- 
sive vertical  stress,  420  and  480  bars,  for  elements  369  and  679,  respect- 
ively, in  Figure  4.11  are  slightly  greater  than  those  at  the  same  points 
from  the  l-MT/300-psi  ( 20. 7-bar) /dense  backfill  calculations  as  shown  in 
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Figures  3.37  and  3.36,  respectively. 


4.4  SHOCK  SPECTRA 

Figure  4.12  presents  the  2 percent  damped  shock  spectra  calculated 
from  the  velocity  wave  forms  in  Figures  4.7  and  4.9.  As  was  observed 
for  the  previous  calculations  (see  Section  3.3),  the  inside  or  the  struc- 
ture roof  (node  513)  has  the  most  severe  calculated  shock  environment, 
and  the  center  of  the  floor  has  the  least  severe  shock  environment  of 
those  points  examined.  Hie  maximum  spectral  response  of  the  structure 
sidewalls  and  roof  occurs  at  frequencies  of  17  to  20  Ha.  The  range  of 
frequencies  for  the  maximum  response  for  the  roof  and  sidewalls  from 
previous  calculations  (see  Section  3.3)  was  from  15  to  17  Hz.  Hie  shock 
spectra  for  all  structural  members  were  made  larger  by  a factor  of  2 or 
less  when  compared  with  those  from  the  10-KT/300-psi  (20.7-bar )/dense 
backfill  case  and  compared  well  with  those  for  the  l-MT/300  psi 
(20.7  bar)/dense  backfill  case. 

4.5  SUMMARY 

The  results  of  the  calculations  presented  in  this  chapter  (i.e., 
nominal  1000-psi  (68.9-bar)  surface  overpressure  loading  from  a 10-KT 
weapon  and  the  dense  backfill  case ) indicated  the  same  trends  as  were 
observed  from  the  calculations  as  presented  in  Chapter  3. 

Hie  maximum  stresses  and  bending  moments  within  and  the  maximum  de- 
flections of  the  structure  for  this  calculation  are  about  the  same  as 
those  from  the  300-psi  (20.7-bar)  1-MT  yield/ dense  backfill  calculation 
presented  in  Chapter  3. 

Hie  results  of  the  shock  spectra  analyses  show  that  the  shock  en- 
vironments for  the  roof,  sidewalls,  and  floor  are  worsened  by  a factor 
of  two  or  less  by  changing  the  nominal  overpressure  from  300  psi 
(20.7  bars)  to  1000  psi  (68.9  bars). 

Hie  predominant  frequency  observed  in  the  motion  time  history  of 
the  center  of  the  roof  was  18  Hz.  It  falls  well  within  the  credible 


frequency  range  as  controlled  by  the  airblast  rise  time  and  duration,  the 
FE  grid,  and  the  controlling  wave  speeds  of  the  backfill  materials.  The 
deflected  shapes  of  the  structure  observed  at  various  times  during  the 
calculation  axe  similar  to  those  observed  in  the  previous  calculations 
carried  out  under  this  project. 


Table  l*.l  Rotations  of  the  Chord  Formed  by  the  'lVo  Bottom  Corners  of 
the  Structure;  10-KT  Yield,  lOOQ-pai  Nominal  Overpressure, 
Dense  Backfill  Calculation . 


■type  Rotation 


Rotation,  rad 


Time , ms 


Counterclockwise  max 
Return  to  upright  position 
Clockwise  max 
Final  (clockwise) 


0.00123 


0.00328 

0.00312 


il 


Table  4.2  Shear,  Throat  and  Bending  Moment  for  Selected  Sections 
Through  the  Structure  at  a Time  of  40  ms  after  the 
start  of  the  10-KT  Dense  Backfill  Calculation. 
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b STRESSES  NORMAL  TO  STRUCTURE 

Figure  4.3  Instantaneous  patterns  of  deflection  of  and 
distribution  of  normal  stresses  on  the  ex- 
terior surfaces  of  the  structure  at  40,  70, 
150,  and  200  ms. 
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Figure  4.5  Moment  and  shear  diagrams  for  tne  structure  sidewalls  at  40  ms  (the  time  that 
the  vertical  stress  at  the  center  of  the  sidewalls  maximized). 
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CHAPTER  5 


PSKUDQSTATIC  CALCULATION 


5.1  CALCULATION  PLAN 


A paeudostatic  calculation  was  conducted  with  the  HONDO  computer 
cod*  for  the  dense  backfill  case,  using  the  same  FE  representation  and 
material  properties  as  discussed  in  Chapter  2.  Hie  calculation  wbb  con- 
ducted with  a linearly  increasing  (with  time)  vertical,  uniformly  dis- 
tributed surface  stress  applied  to  the  top  surface  of  the  FK  representa- 
tion. Hie  Btress  increased  from  0 to  300  psi  (20.7  bars).  A series  of 
one-dimensional  (ID)  calculations  was  conducted  to  determine  the  smallest 
rise  time  and  percentage  of  damping  that  could  be  uBed  and  still  keep 
the  particle  velocities  very  small.  By  trial  and  error,  a rise  time  of 
250  ms  and  a damping  coefficient,  B / * 0.3  , were  determined 

c. 

to  give  acceptable  results.  Hie  pseudostatic  2D  calculation  was  then 
conducted  with  this  rise  time  and  damping  coefficient  and  with  a total 
run  time  of  360  ms.  Hie  300-psi  (20.7-bar)  stress  was  held  constant 
from  250  to  3b0  ms.  Hie  time  step  was  the  same  as  in  previous  problems.** 


5.2  RESULTS 


Except  for  velocity  and  stress  time  histories  for  points  at  and  near 
the  center  of  the  structure  roof,  all  node  point  velocities  became  virtu- 
ally xero  after  350  ms,  while  all  element  stresses  took  on  constant  values 
The  velocity  and  displacement  time  histories  of  node  511 , located  at  the 
center  of  the  structure  roof,  are  shown  in  Figures  5.1  and  5.2,  respec- 
tively. They  present  the  worst  of  the  results  as  far  as  oscillations  in 


The  type  of  damping  used  in  HONDO  and  the  definition  of  the  damping 
coefficients  are  described  in  Reference  2. 

While  this  is  not  a computationally  efficient  way  of  treating  a 
static  problem,  it  turned  out  to  involve  much  less  labor  than  re- 
formulating the  problem  for  a static  code  and,  therefore,  was  more 
economical.  It  had  the  additional  advantage  that  it  ensured  that 
in  dynamic  to  static  comparisons,  "other  things"  were  really  equal. 
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velocity  are  concerned.  In  terms  of  displacements,  the  results  are  quite 
acceptable.  Similarly,  the  horizontal  stress  time  history  shown  in 
Figure  5.3  for  element  511,  located  near  an  extreme  fiber  at  the  center 
of  the  roof,  represents  the  worst  case  for  stress  oscillations  during 
the  pseudostatic  calculation.  It  is  also  an  entirely  adequate  static 
approximation. 

Figure  5.**  shows  the  deflected  shapes  of  structure  and  the  distribu- 
tion of  normal  stresses  on  the  structure  at  80,  160,  and  jbO  ms.  Hie 
vertical  stress  normal  applied  to  the  ground  surface  at  these  times  were 
6.6,  13.2,  and  20.7  bars,  respectively.  However,  the  vertical  stress  in 
the  soil  elements  Just  above  the  central  portion  of  the  structure  at  80 
and  160  ms  was  6.0  and  12.5  bars,  respectively,  while  the  stress  at 
360  ms  was  20.7  bars.  Hie  stress  was  increased  linearly  from  0 to 
20.7  bars  in  250  ms.  Hierefore,  the  stress  was  being  increased  at  80 
and  160  ms  while  the  20.7-bar  stress  had  been  maintained  for  110  ms 
at  the  360-ms  time.  At  the  80-  and  l60-ms  times,  the  differences  in 
applied  stress  at  the  boundary  and  stress  normal  to  the  structure  were 
9.1  and  5.3  percent,  respectively,  Hiese  differences  were  caused  by  the 
fact  that  the  250-ms  rise  time  did  not  allow  the  applied  surface  stress 
to  completely  equalize  throughout  the  problem  due  to  wave  propagation 
effects.  However,  the  differences  (i.e.,  5.3  to  9.1  percent)  are  not 
considered  significant. 

As  can  be  seen  in  Figure  5.**»  the  stress  and  deflection  patterns  are 
symmetrical  about  a vertical  centerline  through  the  structure.  Hie  normal 
stress  on  the  roof  of  the  structure  is  uniform  at  each  stress  level  except 
over  the  sidewalls  where  the  stress  is  increased  due  to  negative  arching. 
Hie  stresses  on  the  sidewalls  increase  with  increase  in  applied  surface 
stress.  Hie  stress  distributions  normal  to  the  sidewall  are  concave  to 
the  structure  sidewall  and  maximize  near  the  center  of  the  sidewalls. 

A portion  of  the  stress  on  the  sidewalls  is  due  to  partial  passive  earth 
pressure  caused  by  the  sidewalls  bowing  out  against  the  backfill.  Hie 
stresses  normal  to  the  base  of  the  structure  show  a concentration  of 
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stress  under  the  sidewalls  and  are  reasonable  when  compared  with  the 
deflected  shape  of  the  structure.  The  patterns  of  deflected  shape  of 
and  the  stresses  normal  to  the  structure  at  3b0  ms  are  similar  to 
those  from  previous  calculations  at  the  time  that  the  deflection  of 
the  center  of  the  roof  maximized  as  shown  in  Figures  3.6  and  4.3,  The 
maximum  downward  deflection  of  the  center  of  the  roof  for  the  pseudo- 
static calculation  is  10.5  centimetres,  and  the  maximum  relative  deflec- 
tion of  the  center  of  the  roof  (i.e.,  the  deflection  of  the  center  of 
the  roof  with  respect  to  a chord  drawn  through  the  top  corners  of  the 
structure)  is  4.4  centimetres. 

Axial  thrust,  3hear  force,  and  bending  moment  per  unit  width  of 
the  atructure  were  calculated  at  a time  of  360  ms  and  are  tabulated  in 
Table  5.1  for  the  28  structural  sections  shown  in  Figure  2.20.  Moment 
and  shear  diagrams  are  plotted  for  the  structure  roof  and  floor  in 
Figure  5.5  and  for  the  structure  sidewalls  in  Figure  5.0.  The  moments 
in  the  roof  and  floor  appear  to  maximize  near  the  center  and  ends  of 
the  spans  indicating  some  degree  of  fixity  due  to  the  sidewalls. 

Examination  of  Figures  5.1  through  5.6  indicates  that  results  ob- 
tained from  the  pseudostatic  HONDO  calculation  are  reasonable  approxi- 
mations to  a true  static  calculation.  The  results  of  this  calculation 
will  be  used  in  the  next  chapter  to  establish  ratios  of  dynamic  to 
static  response. 
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Table  5.1  Shear,  Thrust, and  Bending  Moment  for  Selected  Sections 
Through  the  Structure  at  a Time  of  360  ms  after  the 
start  of  the  10-KT  "Pseudo"  Static  Dense  Backfill 


Calculation . 
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SHEAR  FORCE 

AXIAL  THRUST 

9EN0ING  MOMENT 
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AA 
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5.3  Horizontal  stress  tixae  histories  for  element  511  fron  the  oseudostatic  calculation 
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Figure  5»**  Normal  stress  and  deformation  pattern  for 
three  applied  stress  levels  from  the 
pseudostatic  calculation. 
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Figure  5.5  Moment  and  shear  diagrams  for  the  structure  roof 
and  floor  at  360  ms  from  the  pseudostatic  calcu- 
lation. 


CHAPTER  b 


ANALYSIS 


INTRODUCTION 


As  discussed  in  Chapter  3 a change  in  backfill  properties  from  those 
of  a poorly  compacted  material  to  those  of  the  same  material  with  a high- 
er degree  of  compaction  reduces  the  calculated  load  on  and  response  of  a 
shallow-buried  box  structure  subjected  to  surface  airblast  loading.  Idle 
first  part  of  this  chapter  explores  why  this  is  so.  The  second  portion 
of  this  chapter  deals  with  the  effects  of  overpressure  (P  ) and  yield 

(W)  variations.  Obviously,  as  P and  W increase,  so  do  the  maxi- 

so 

mum  amplitudes  of  structural  response.  'Ihis  is  well  documented  from 
linear  elastic  dynamic  structural  analyses  of  protective  structures. 

The  relatively  small  increase  in  response  due  to  a 3.33  increase  in  P 

SO 

indicated  in  this  study  was  considered  surprising,  and  an  effort  was  made 
to  try  to  understand  the  phenomena  that  were  responsible  because,  if  true 
this  trend  has  implications  on  the  accuracy  required  in  targeting  with 
low-yield  weapons. 


6.2  PHENOMENA  RESPONSIBLE  FOR  DIFFERENCES  CAUSED  BY  A CHANGE  IN  BACKFILL 
PROPERTIES 


Ihe  calculations  described  in  Chapter  3 were  conducted  for  two  back' 
fill  conditions  (loose  and  dense)  and  for  two  overpressure  loading  rela- 
tions (i.e.,  nominal  300-psi  (20.7-*bar)  airblast  overpressure  loading 
from  1-MT  and  10-KT  weapons).  As  far  as  response  of  the  linear  elastic 
idealized  structure  is  concerned,  the  conclusion  that  may  be  drawn  is 
that  the  effect  of  using  loose  rather  than  dense  backfill  properties  for 
either  weapon  yield  is  to  increase  the  calculated  deflections  of  and  the 
stresses  within  the  structure  by  25  to  Lo  percent.  The  purpose  of  this 
section  is  to  examine  the  phenomena  that  could  be  responsible  for  these 
differences. 

Tdie  deflections  of  the  center  of  the  roof,  the  moments  in  the  roof. 


and  the  bowing  out  of  the  sidewalls  of  the  structure  were  greater  for 
the  loose  backfill  conditions.  Hie  dense  backfill  appears  to  cause  more 
restraint  on  the  structure  sidewalls;  therefore,  the  structure  roof  be- 
haves more  like  a fixed  beam  when  surrounded  by  dense  backfill  and  more 
like  a simply  supported  beam  when  surrounded  by  loose  backfill.  Hie 
greater  lateral  restraint  on  the  sidewalls  for  dense  backfill  conditions 
could  be  caused  by  higher  radial  stresses  transmitted  to  the  structure 
through  the  dense  backfill  than  through  the  loose  backfill  and/or  by 
its  greater  strength  and  stiffness,  which  would  provide  larger  partial- 
passive  earth  pressure  per  unit  outward  wall  deflection. 

6.2.1  Vertical  Stress  Attenuation 

The  comparison  of  stress-strain  properties  for  a state  of  UX  for 
the  loose  and  dense  backfills  shown  in  Figure  2,l4  indicates  that  much 
more  energy  will  be  dissipated  as  the  applied  stress  pulse  passes  through 
the  loose  backfill  than  through  the  dense.  Hiis  is  due  to  greater  hys- 
teresis in  the  loose  backfill  stress-strain  relation.  Hie  practical 
impact  of  this  is  that  the  vertical  airblast -induced  stress  should  be 
attenuated  more  as  it  passed  through  the  loose  backfill  material  than 
through  the  dense  backfill. 

Figure  6.1  shows  the  variation  of  peak  vertical  stress  with  depth 
in  the  backfill  from  the  1-MT  calculations  at  a point  3.5  metres  up- 
stream of  the  leading  edge  of  the  structure.  As  also  noted  in  the 
figure,  the  maximum  vertical  stresses  are  roughly  the  same  in  the  two 
calculations, which  was  not  expected.  Hie  discrepancy  may  be  explained 
in  the  following  manner.  Since  the  airblast  positive  phase  durations 
were  sufficiently  long  with  respect  to  the  depth  of  cover  and  the  wave 
propagation  properties  of  both  backfill  materials,  the  ID  wave  propagation 
theory  for  linear  hysteretic  materials  (Figure  6.2)  indicates  that  stress 
attenuation  at  the  depth  of  the  roof  should  have  been  negligible  in  both 
cases.  Hiis  is  in  agreement  with  the  results  presented  in  Figure  6.1. 
Hence,  both  theory  and  calculation  results  indicate  that  the  maximum 
vertical  free-field  stress  at  the  depth  of  the  roof  was  almost  the  same 
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for  the  two  backfill  materials.  Since  the  constitutive  model  is  isotropic, 
it  may  be  concluded  that  radial  stress  attenuation  is  not  the  primary 
cause  of  reduced  lateral  restraint  on  the  sidewalls  for  the  loose  back- 
fill condition. 

6.2.2  Lateral  at  Rest  Earth  Pressures 

In  addition  to  dixlferences  in  stress-strain  behavior,  the  backfills 
differ  in  their  UX  stress  path  behavior  (Figures  2.11  and  2. IS);  however, 
this  difference  is  small.  Had  the  differences  in  the  loading  stress  paths 

(whose  slopes  are  ~ ) been  significant,  the  at-rest  lateral  pres- 

sure  coefficients  K = — ^r-lu?.n,'a.  ■ = - — - — would  have  been  different, 

o O ...  1 - V 

vertical 

and  equal  vertical  stresses  adjacent  to  the  structure  would  have  produced 
different  lateral  stresses  on  a rigid  (nonmoving)  structure.  For  the 

loose  and  dense  backfills  over  the  0-  to  20-bar  range,  the  K 's  are  0.50 

o 

and  0.4U,  respectively.  This  difference  seems  too  small  to  have  any 
significant  effect  on  structure  loading. 

6.2.3  Lateral  Passive  Earth  Pressures 

Figure  6.3  shows  horizontal  stress  time  histories  for  elements  ad- 
jacent to  the  blastward  structural  sidewall  from  the  l-MT/300-psi 
(20.7-bar )/loose  and  dense  backfill  calculations.  As  also  noted  in  the 
figure,  the  horizontal  stress  during  the  50-  to  75-ms  time  range  is  much 
greater  for  the  dense  backfill  case  than  for  the  loose  backfill  case 
in  spite  of  the  fact  that  K0's  are  similar.  This  time  range  corresponds 
with  the  times  at  which  the  maximum  downward  deflection  of  the  roof  and 
the  maximum  outward  deflection  of  the  sidewalls  occur  as  shown  in  Fig- 
ure 3.6.  This  also  is  the  time  that  the  maximum  moments,  thrusts,  and 
shears  occur  in  the  structure  as  shown  in  Figures  3.11  through  3.13. 

Figure  6.3  is  consistent  with  the  idea  that  the  structure  is  being  re- 
strained more  in  the  horizontal  direction  for  the  dense  backfill  case. 

This  restraint  would  tend  to  cause  the  roof  and  floor  to  have  a greater 
degree  of  end  fixity.  This  would,  in  turn,  account  for  the  fact  that 
the  deflections  of  and  stresses  within  the  structure  roof  and  floor  are 


less  for  the  dense  backfill. 

After  examining  the  wave  speeds  for  the  backfill  material  shown  in 
Table  6.1,  it  is  apparent  that  the  stress  is  arriving  for  both  backfill 
cases  faster  than  it  is  possible  for  the  stress  to  propagate  downward 
through  the  backfill  materials.*  After  examining  various  paths  for 
first  arrivals, the  conclusion  was  drawn  that  the  path  for  the  fastest 
arrival  of  stress  at  the  raiddepth  of  the  backfill  adjacent  to  the  struc- 
ture is  downward  through  the  structure  sidewall.  The  time  for  the  stress 
to  arrive  at  the  center  of  element  276,  located  at  a depth  of  6.5  metres 
below  the  ground  surface  and  0.3Q  metres  from  the  structure  sidewall  by- 
taking  a path  through  the  structure  sidewall,  is  shown  in  Table  6.2  for 
each  of  the  backfill  cases.  The  initial  loading  wave  speeds  were  used 
for  the  backfill  materials  and  the  concrete.  The  arrival  times  at 
element  278  computed,  as  shown  in  Table  6.2,  for  the  J00  psi/l-MT  weapon 
yield/loose  and  dense  backfills  are  1*5.9  and  37.5  ms,  respectively. 

These  times  compare  favorably  with  the  code  calculated  first  arrival 
times  of  1*8.0  and  39.0  ms  shown  in  Figure  6.3.  Lines  showing  the  cal- 
culated times  of  arrival  of  stress  through  the  backfill  are  presented 
in  Figure  6.3.  They  show  that  the  stress  propagating  vertically  down- 
ward through  the  backfill  is  arriving  much  later  than  that  which  is 
coming  through  the  structure.  Figure  6.1*  shows  the  same  trends  for 
the  10-KT  loading  case.  From  these  analyses,  it  was  concluded  that  much 
of  the  horizontal  stress  experienced  by  the  elements  shown  in  Figure  6.^ 
was  due  to  causes  other  than  wave  propagation  through  the  backfill.  As 
the  sidewall  bows  out  againBt  the  backfill  due  to  loading  on  the  struc- 
ture roof,  passive  earth  pressure  is  developed.  The  passive  resistance 
is  much  greater  for  the  dense  backfill  case  during  the  time  that  the  maxi- 
mum deflections  and  moments  in  structural  sections  occur  (Figure  6.3) 
and  accounts  for  the  greater  structural  deflections  and  moments  for  the 
loose  backfill  cases. 


The  code  and  model  are  not  in  error;  this  was  proven  by  the  ID  cal- 
culations in  which  the  expected  propagation  velocities  were  observed. 
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Figure  b.5  shows  the  shear  failure  surfaces  for  the  two  backfill  mu- 
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terials.  At  a given  mean  normal  sti'ess  level,  the  dense  backfill  has 
roughly  twice  the  strength  of  the  loose  backfill.  Since  the  muximum  pos- 
sible passive  earth  pressure  that  can  be  mobilized  in  a given  material, 
increases  with  strength,  the  trends  in  stresses  in  Figures  6.j  and  b.4 
are  consistent  with  passive  earth  pressure  phenomena. 

6.2.4  Stress  Wave  Reflections  off  the  Concrete 

Another  phenomenon  that  causes  differences  in  structural  response 
for  the  loose  and  dense  backfill  cases  is  reflection  of  stress  on  the 
roof  due  to  the  variation  in  stress-strain  properties  (Figure  2.14)  of 
the  backfill  materials.  Figure  b.b  shows  the  vertical  stress  time 
histories  for  points  Just  above  the  structure  roof  in  the  loose  and  dense 
backfill  calculations  for  1-MT  and  10-KT  weapons.  Also,  in  Figure  b.b, 
the  vertical  stress  for  both  weapon  yields  and  backfill  conditions ' is 
highly  oscillatory.  This  phenomenon  is  caused  by  reflection  off  the 
relatively  stiff,  thick  roof  structure  at  a depth  of  1.52  metres  and 
by  the  stress  wave  reflection  at  the  free  surface,  i.e.,  the  ground 
surface.  The  period  of  oscillation  should  be  approximately  411/C^  » 
where  H is  depth  to  the  top  of  the  structure  and  C is  the  unloading- 
reloading wave  velocity  (Table  6.1 ).  For  the  loose  and  dense  backfill 
curves,  the  periods  calculated  by  the  formula  are  10  and  7 ms,  respec- 
tively. The  average  periods  in  Figure  b.b  are  9-5  and  8,5  ms  for  the 
loose  and  dense  backfills,  respectively.  This  oscil-lary  stress  time 
history  phenomenon  has  been  seen  experimentally  in  the  WES  small  blast, 
load  generator  (Reference  10).  Since  there  is  a greater  impedance  mis- 
match between  the  loose  backfill  and  the  concrete  than  between  the  dense 
backfill  and  the  concrete,  the  reflected  vertical  stresses  for  the  loose 
backfill  cases  are  greater  than  that  for  the  dense  backfill  cases  (Fig- 
ure 0.6).  Hence,  for  the  same  applied  airblast  loading  at  the  ground 
surface,  a greater  vertical  loading  is  experienced  by  the  structure  roof 
for  the  loose  backfill  cases  than  for  the  respective  dense  backfill  cases. 
It  is  concluded  that  these  differences  in  reflection  phenomena  coupled 
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with  differences  in  passive  earth  pressure  effects  are  the  primary  causes 
of  larger  structural  deflections  and  moments,  thrusts,  and  shears  within 
the  structures  for  the  loose  backfill  cases. 

6.3  EFFECT  OF  SOIL  COVER  ON  NATURAL  FREQUENCIES 

In  Appendix  A,  eigenvalue  analyses  of  the  box  structure  both  in  vacuo 
and  embedded  in  elastic  media  are  reported.  Figure  6.7  shows  the  lower 
two  calculated  frequencies  and  corresponding  mode  shapes  for  the  structure 
calculated  during  the  eigenvalue  analyses  for  the  in  vacuo  and  embedded 
(stiff)  conditions.  Summarized  below  are  the  results  for  the  lower  modes: 


Case  Studied  Hz  Hz  Remarks 


In  Vacuo  15.1  18.2 

Itabedded  (Stiff)  12.0  16.7  Based  on  unloading  moduli  for  earth 

media 

kknbedded  (Soft)  5.1  0.8  Media  10  times  softer  than  in  stiff 

media  embedded  analysis 

(Added  Mass]  (11. 7]  [14.2] 

ilnbedment  causes  a decrease  in  frequency;  the  softer  the  media  the  larger 
is  the  decrease.  In  simple  terms,  this  decrease  is  occurring  because  the 
media  contributes  more  to  the  mas3  of  the  system  than  to  its  stiffness. 
Reference  5 describes  a way  of  estimating  the  added  mass  effect.  Its  ap- 
plication to  the  in  vacuo  modes  is  shown  at  the  bottom  of  the  tabulation 
above, and  it  produces  results,  which  are  in  fairly  good  agreement,  with 
the  embedded  analysis  that  employed  elastic  constants  based  on  the  un- 
loading properties  of  the  media,  'llie  deformed  shapes  of  the  structure 
(Figure  3.6)  at  times  of  maximum  response  appear  to  be  a combination  of 
the  side  sway  and  the  roof-floor  bending  modes  with  the  latter  being  pre- 
dominant. 

6.4  APPARENT  DAMPING  IN  EQUIVALENT  S INGLE-DBGRKE-OF-FREKDOM  SYSTEM 

One  elementary  way  of  considering  the  dynamic  response  of  the  linear 


Bide 

Sway 

Mode 

Hz 

Roof-Floor 

Bending 

Mode 

Hz 

Remarks 

15.1 

18.2 

12.0 

16.7 

Based 

media 

on  unloading  moduli  for 

earth 

5.1 

b.8 

Media 

media 

10  times  softer  than  in 
embedded  analysis 

stiff 

(11.71 

[14.2] 
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elastic  structure  studied,  in  this  investigation  is  by  means  of  a linear 
elastic  single-degree-of-freedom  (BDF)  system  analysis.  If  a single  mode 
of  the  2D  structure  is  predominant  in  the  2D  response  calculations  (and 
this  is  the  case),  an  equivalent  SDF  model  can  be  created,  which  has  the 
natural  frequency  of  this  mode.  Figure  6.8  shows  the  roof  centerline 
deflection  time  histories  for  four  of  the  five*  dynamic  2D  calculations 
and  the  maximum  deflection  in  the  300-psi  (20.7-bar),  pseudostatic  sur- 
face load  calculation.  Also  presented  in  the  figure  are  the  surface 
overpressure  loading  functions  used  in  these  calculations.  From  the 
deflection  time  histories,  the  periods  of  the  oscillation  can  be  scaled: 

Time  from 

Airblast  Period  of 
Arrival  First 

Yield  to  Maximum  Post-Peak. 

W Displacement  Oscillation 

MT  Backfill  ms ms 

1 Dense  44  (ll  Hz)  52  (19  Hz) 

1 Loose  46  (ll  Hz)  65  (l6  Hz) 

0.01  Dense  40  (13  Hz)  54  (19  Hz) 

0.01  Dense  32  (lb  Hz)  52  (19  Hz) 

The  post-peak  frequency  in  the  dense  backfill  cases  are  in  fairly  good 

agreement  with  the  embedded  (stiff)  and  in  vacuo  roof-floor  bending  modes 

In  the  loose  backfill  case,  the  frequency  is  lower  than  in  the  dense,  but 

the  percentage  change  is  much  smaller  than  in  the  eigenvalue  analyses 

discussed  in  Section  0.3. 

In  Figure  6.9,  the  deflection  of  the  centerline  of  the  roof  relative 
to  the  average  of  the  corners  of  the  floors  is  plotted  for  the  1-MT/ 
300-psi  (20.7-bar )/dense  backfill  calculation.  When  the  deflection  of 
the  underlying  media  is  subtracted  out,  the  relative  deflection  curve  is 
seem  to  be  oscillating  about  a curve  similar  in  shape  to  the  applied  air- 
blast  loading  function  with  a period  slightly  longer  than  that  observed 


Peak 

Overpressure 


Ihe  10-KT  loose  backfill  calculation  was  not  considered  here  because 
cursory  examination  showed  it  would  show  similar  loose  and  dense  trends 
as  the  1-MT  calculations. 


[ 
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in  the  absolute  deflection  curve.  All  of  the  roof  centerline  deflection 
curves  in  Figure  6.8  were  processed  in  this  manner  to  obtain  the  periods 
given  in  Table  6.3. 

The  periods  for  the  dense  backfill  cases  in  Table  6.3  are  approxi- 
mately 60  ms.  A SDK  system  vith  an  undamped  natural  period  of  60  ms 
was  subjected  to  a unit  load  whose  time  variations  were  the  same  as  the 
overpressure  functions  shown  in  Figure  6.8.  Dynamic  responses  were  calcu- 
lated for  the  cases  of  0,  20,  and  40  percent  of  critical  damping.  The 
dynamic  load  factors  (DLF)  from  these  calculations  (i.e.,  the  ratio  of 
maximum  deflection  of  the  SDF  system  to  the  deflection  of  the  same 
system  under  the  same  peak  load  applied  statically)  are  given  in 
Table  6.3  and  can  be  seen  to  decrease  significantly  as  the  percent  of 
critical  damping  increases.  These  DLF  were  compared  with  2D  DLF  derived 
from  the  HONDO  code  calculations  by  the  following  formula: 

maximum  deflection  of  centerline  of  the  roof  relative  to  the  floor 
corners  in  dynamic  HONDO  calculation 


2D  DLF  = 


relative  deflection  at  the  same  surface  stress  level  from  the 
pseudostatic  calculation 


Since  the  structure  itself  is  linear  elastic,  this  comparison  is  legitimate 
and  shows  that  in  order  to  achieve  responses  similar  to  the  2D  calculations, 
the  SDF  system  must  be  20  to  40  percent  critically  damped  with  higher  per- 
centage of  damping  being  appropriate  for  the  shorter  duration  (lower  yield) 
external  loadings. 

Since  there  i3  no  source  of  damping  within  the  structure,  what  is  the 
mechanism  present  in  the  2D  calculations  that  has  the  appearance  of  equiv- 
alent viscous  damping?  There  are  two  possible  sources:  the  energy  ab- 
sorbing (hysteresis)  capacity  of  the  surrounding  inelastic  media  and 
radiation  "damping,"  i.e.,  the  propagation  of  energy  away  from  the  loaded 
boundary  due  to  the  presence  of  a large  linear  elastic  (or  inelastic) 
region.  Reference  11  concludes  that  "all  footing  soil  systems  are  strongly 
damped"  because  of  the  propagation  of  energy  into  the  subsoil  even  if  the 
subsoil  behaves  as  a linear  elastic  medium.  Ihe  box  structure  here  can 
be  viewed  as  a large  mat  footing  sitting  on  a foundation  material,  and 
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its  "rigid  body"  motion  should  be  expected  to  be  strongly  damped.  In 
fact,  it  is  as  shown  by  the  dotted  curve  in  Figure  6.9.  Hie  same  phe- 
nomenon should  have  an  effect  on  the  nonrigid  body  modes  of  the  structure 
and  is  certainly  responsible  for  some  of  the  apparent  damping  in  the 
roof— floor  bending  mode  seen  in  the  calculated  structural  response 
(References  12  and  13).  Hie  hysteresis  of  the  backfill  should  also 
cause  absorption  of  energy.  Neither  of  these  factors  are  usually  con- 
sidered in  SDF  system  analyses  of  shallow-buried  protective  structures 
subjected  to  airblast  loadings. 

While  this  study  is  not  conclusive,  it  strongly  suggests  that  SDF 
system  models  used  in  targeting  analyses  of  shallow-buried  structures 
should  include  a dissipative  mechanism  to  account  for  the  energy  re- 
moved from  the  system  via  radiation  and  hysteresis.  Hie  appropriate 
damping  appears  to  exceed  20  percent  of  critical  for  the  class  of  prob- 
lems studied  here.  In  cases  where  the  overpressure  decays  very  slowly 
(the  combination  of  low  ( 100  psi)  surface  overpressures  and  high 
yields),  ignoring  this  mechanism  has  little  effect  on  the  maximum  re- 
sponse of  the  system.  For  high  overpressure,  and/or  low  yields,  the 
decay  time  to  one-half  the  peak  overpressure  can  be  less  than  the 
natural  period  of  the  system,  and  the  system  response  can  be  strongly 
affected  by  the  presence  of  dissipative  mechanisms.  Hie  presence  of 
these  dissipative  mechanisms  is  thought  to  be  the  reason  why  the 
differences  in  response  of  the  2D  elastic  structure  due  to  a factor 
of  3 variation  in  overpressure  were  not  as  large  as  expected  based  on 
previous  thinking  (largely  conditioned  by  undamped  system  analyses) 
about  the  problem. 

As  an  example  of  the  potential  implications  consider  the  following 
tabulation : 


Overpressure 
psi  (bar) 

Yield 

MT 

Distance 

from 

GZ 

km 

Maximum 
Relative* 
Deflection 
of  Roof 
Based  on 

0*  of  Damped 
SDF  Model 

cm 

Actual 
HONDO  Code 
Relative 
Deflection 

cm 

Overestimate 
of  Structural 
Deformation  by 
Undamped 

SDF  System 
% 

300  (20.7) 

.01 

.15 

5.3 

3.4 

71 

1000  (68.9) 

.01 

.09 

10.9 

5.6 

95 

300  (20.7) 

1.00 

.67 

6.1 

5.3 

15 

Die  undamped  SDF  system  for  the  two  10-KT  cases  overestimated  the  struc- 
ture roof  deflection  by  a factor  of  nearly  2.  If  6 centimetres  of  de- 
flection were  assumed  to  be  the  amount  required  to  cause  failure ,##  one 
would  have  to  achieve  a 1000-psi  (68.9-bar)  overpressure  over  the  struc- 
ture to  cause  the  failure  that  the  SDF  analysis  would  predict  for  the 
300  psi  (20.7  bars).  Hie  CiSP  would  have  to  be  reduced  by  a factor  of 
1.7  to  maintain  the  same  single-shot  kill  probability  for  the  10-KT 
yield.  By  introducing  the  appropriate  percentage  of  damping  into  the 
analysis,  this  trend  becomes  apparent. 


* DLF  of  0 percent  damping  x P /300  psi  x static  deflection  at  300  psi. 

SO 

**  For  a real  structure,  this  would  be  a ridiculously  small  value. 
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Figure  6.1  Maximum  vertical  stress  versus  depth 
for  a station  3.5  m upstream  of  the 
structure,  1-MT  case. 
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Stress  attenuation  in  linear  hysteretic  media 
similar  to  loose  and  dense  backfills  for 
20-bar  airblast  loadings . 
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Figure  6.7  Side  away  and  roof-floor 
bending  modes. 
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Figure  6.9  Relative  roof  deflection  versus  time  for  300-psi  ( 20-bar )/l-MT/ dense 
backfill  calculation. 


CHAPTER  7 


SUMMARY,  CONCLUSIONS,  AND  RECOMMENDATIONS 


7.1  SUMMARY 


Five  dynamic  and  one  pseudostatic  plane  strain,  FE  calculations  were 
performed  with  the  HONDO  code.  In  the  dynamic  calculations,  the  response 
of  a shallow-buried  linear  elastic  box  structure  was  examined  for  several 
conditions.  Nonlinear  backfill  properties  were  varied  from  those  of  a 
dense,  well-compacted  to  those  of  a loose,  poorly  compacted  material. 

Both  the  backfill  and  underlying  materials  were  modeled  as  nonlinear 
and  repetitive  using  a cap-type  constitutive  model.  Peak  overpressure 
and  yield  associated  with  the  traveling  surface  airblast  load  were  also 
varied. 


7.2  CONCLUSIONS 


For  the  idealized  2D  plane-strain  problem  studied,  the  following 
conclusions  can  be  drawn: 

1.  As  the  stress-strain  and  strength  properties  of  the  backfill 
material  are  changed  to  those  of  a loose  shale  backfill  from  those  of 
the  same  material  in  a well-compacted  state,  the  loads  on,  deflections 
of,  and  shears,  moments,  and  thrusts  within  the  2D  elastic  structure 
increase  by  25  to  40  percent. 

2.  The  primary  causes  of  this  increase  are  (1)  a decrease  in 
corner  fixity  resulting  from  decreased  passive  earth  pressure  effects 
along  the  sidewalls,  and  (2)  higher  reflected  stresses  at  the  roof- 
backfill  interface  due  to  the  greater  impedance  mismatch  in  the  loose 


backfill  case. 

3.  As  yield  is  increased  from  10  KT  to  1 MT,  an  idealized  structure 
at  the  300-psi  (20.7-bar)  overpressure  level  undergoes  an  increase  in 
peak  response  (in  terras  of  maximum  applied  stress,  deflections,  and 
internal  moments,  thrusts,  and  shears)  of  approximately  45  to  60  percent. 

4.  An  increase  in  peak  overpressure  from  300  to  1000  psi  (20.7  to 
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68.9  bars)  at  a constant  10-KT  yield  generally  resulted  in  a 1*0  to  60 
percent  increase  in  structural  response. 

5.  Supporting  linear  elastic  eigenvalue  analyses  revealed  mode 
shapes  and  frequencies  consistent  with  those  observed  in  the  HONDO 
calculations.  Hie  decrease  in  backfill  stiffness  caused  a noticeable 
decrease  in  frequency  associated  with  the  predominant  mode  of  deformation. 

6.  The  linear  SDF  system  that  best  simulated  the  response  of  the 
2D  linear  elastic  structure  to  the  locally  applied  airblasts  was  one 
whose  natural  frequency  was  the  same  as  that  of  the  predominant  mode  of 
deformation  and  that  was  damped  at  20  to  40  percent  of  critical  damping. 
Hie  latter  appears  necessary  to  acccrnt  for  energy  lost  to  the  surrounding 
media. 

7.3  RECOMMENDATIONS 

In  targeting  analyses  of  shallow-buried  box  structures,  it  is  con- 
ventional to  use  an  equivalent  SDF  model  of  the  facility  (References  5 
and  9 provide  examples  of  SDF  models)  and  to  assume  that  the  roof  load- 
ing is  the  surface  overpressure  produced  by  a given  detonation.  The 
studies  performed  here  suggest  that  this  may  result  in  an  underestimation 
of  the  hardness  of  the  facility  if  the  airblast  has  a short  effective 
duration  because  the  radiation  damping  provided  by  the  surrounding  media 
and  its  hysteresis  both  permit  energy  to  be  removed  from  the  structure. 
Apparent  damping  of  more  than  20  percent  of  critical  was  noted  in  the 
cases  studied  here  and  was  responsible  for  reduced  deformation  levels 
under  short-duration  transient  loading.  In  view  of  the  potential  im- 
plications of  this  on  the  delivery  accuracy-yield  trade-off,  further 
studies  of  apparent  damping  should  be  undertaken  to  determine  the  class 
of  cases  for  which  it  is  significant  and  how  best  to  account  for  the 
additional  hardness  provided  by  the  surrounding  media's  dispersive  and 
dissipative  nature. 


- - 
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APPENDIX  A 


IN  VACUO  AND  EMBEDDED  MODAL  ANALYSES 
FOR  A HYPOTHETICAL  STRUCTURE  IN  PLANE  STRAIN 

A.l  OBJECTIVE  AND  SCOPE 

Prior  to  performing  FE  calculations,  a series  of  in  vacuo  modal 

II- 

analyses  was  conducted  with  the  SAP  IV  elastic  FE  code  to  (l)  deter- 
mine the  effects  on  the  FE-calc ulated  natural  frequencies  of  the  hy- 
pothetical structure  using  constant  strain  elements  (i.e.,  the  type 
employed  by  the  HONDO  code)  instead  of  higher  order  elements,  and  (2) 
aid  in  determining  the  minimum  number  of  elements  needed  across  struc- 
tural sections  in  order  to  adequately  model  the  response  of  the  struc- 
ture. A series  of  embedded  modal  analyses  was  subsequently  performed 
to  (l)  determine  the  effect  of  material  property  variations  on  the 
natural  frequencies  of  the  buried  structure,  and  (2)  assist  in  planning 
the  follow-on  structure/backfill  interaction  calculations  so  as  to  have 
a sufficient  range  of  frequency  response  to  permit  an  adequate  assess- 
ment of  the  actual  stress  and  motion  patterns  under  the  postulated  air- 
blast  environments. 

The  purpose  of  this  appendix  is  to  present  the  results  of  both  the 
in  vacuo  and  the  embedded  modal  analyses  and  to  discuss  their  implication 
with  regard  to  the  overall  study. 

A. 2 IN  VACUO  ANALYSES 

Since  SAP  IV  cannot  rigorously  handle  in  vacuo  analyses,  an  approxi- 
mation was  made  by  placing  very  soft  springs  at  the  corners  of  the  struc- 
ture in  such  a way  that  the  rigid  body  modes  (i.e.,  those  that  are  a 
function  of  the  boundary  springs  only)  did  not  interfere  with  the  pure 
in  vacuo  or  deformation  modes  of  the  structure.  The  structure  and  the 
springs  attached  to  it  are  shown  in  Figure  A.l;  the  concrete  was  assumed 
to  have  a Young's  modulus  E « 2 x 10  psi  (1.38  x 10J  bars)  and  a 
Poisson's  ratio  v = 0.25.  Eigenvalue  calculations  for  a two-element- 
per-section  structural  idealization  were  conducted  using  10-pounds /inch 
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(1751.2-Newton/metres)  and  100-pounds /inch  (17,512-Newton/metres)  spring 
constants.  Of  course,  the  three  rigid  body  frequencies  were  significantly 
altered  by  the  spring  constant  variation,  but  the  frequencies  computed 
for  the  first  ten  structural  or  deformation  modes  were  identical  to  four 
significant  figures.  Thus,  the  spring  constants  for  subsequent  in 
vacuo  analyses  were  all  100  pounds/inch  (17»512  Newton/metres) . 

Additional  calculations  were  then  performed  for  structure  idealiza- 
tions employing  four  and  eight  elements  per  section;  these  FE  descrip- 
tions of  the  structure  are  shown  along  with  the  previously  used  two- 
element-thickness  description  in  Figure  A. 2.  Ihe  elements  used  in  all 

of  the  calculations  described  thus  far  were  of  the  higher  order  type 

1U 

developed  by  Wilson.  These  higher  order  elements  are  provided  in 
SAP  IV  by  adding  incompatible  displacement  modes  to  the  basic  quadri- 
lateral element  description  in  order  to  improve  the  bending  properties 
of  the  element.  Hie  frequencies  calculated  for  the  first  ten  struc- 
tural modes  are  given  in  Table  A.l;  those  for  the  eight-element  descrip- 
tion are  assumed  to  best  approximate  the  true  frequencies  for  this 
structure.  By  examining  Table  A.l,  it  was  found  that  the  frequencies 
for  the  two-  and  four-element-thickness  idealizations  differ  only 
slightly  from  those  computed  for  the  more  detailed  eight-element 
idealization. 

The  HONDO  code  being  used  for  the  structure/backfill  interaction 
calculations  only  employs  isoparametric  quadrilateral  elements  and  as- 
sumes that  the  strains  computed  at  the  center  of  the  element  represent 

p 3 

average  values  for  the  entire  element.  * Such  elements  are  hereinafter 
referred  to  as  constant  strain  elements  and  were  used  in  another  series 
of  two-,  four-,  and  eight-element-sectional-thickness  calculations. 

The  frequencies  obtained  for  the  first  ten  structural  modes  for  the 
constant  strain  element  in  vacuo  calculations  are  given  in  Ihble  A. 2. 

As  with  the  higher  order  elements,  there  is  little  difference  in  the 
two-,  four-,  and  eight-element-thickness  frequencies;  and  the  structural 
mode  shapes  were  similar  for  all  six  of  these  calculations.  However,  a 
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comparison  of  the  frequencies  in  Table  A. 2 with  those  in  Table  A.l  shows 
that  those  obtained  using  the  hi$ier  order  element  description  are  not 
significantly  improved. 


The  four-element-thickness  structural  approximation  was  chosen  for 
both  the  structure/backfill  interaction  parameter  study  and  the  em- 
bedded eigenvalue  analyses,  because  the  in  vacuo  results  indicate  that 
the  rather  large  additional  costs  associated  with  using  eight-element 
thicknesses  are  not  justifiable.  Since  the  structure/backfill  interaction 
calculations  will  have  to  employ  constant  strain  elements,  such  elements 
were  selected  for  the  embedded  eigenvalue  analyses.  Thus,  for  purposes 
of  subsequent  comparison,  the  frequencies  and  deflected  structure  shapes 
for  the  three  rigid  body  modes  and  the  first  10  structural  modes  obtained 
from  the  in  vacuo  analysis  using  a constant  strain  element  description 
and  the  four-element-thickness  approximation  are  given  in  Figure  A. 3. 


A. 3 EMBEDDED  ANALYSIS 

Two  embedded  eigenvalue  analyses  were  conducted,  one  using  a lower 
bound  approximation  to  the  stiffness  of  the  soils  surrounding  the  struc- 
ture and  the  other  an  upper  bound  approximation.  The  FE  grid  represen- 
tation for  the  plane-strain  structure/backfill  interaction  calculation 
study  is  shown  in  Figure  A. 4 along  with  the  portion  of  grid  included 
in  the  embedded  modal  analyses  (i.e,,  area  ABCD).  The  boundary  constraints 
for  the  SAP  IV  modal  analyses  were  similar  to  those  planned  for  the  HONDO 
study,  i.e.,  boundaries  AB,  BC,  and  CD  were  rollers,  boundary  AD  was 
free,  and  points  B and  C were  fixed.  The  free-field  material  descrip- 
tions as  well  as  those  for  the  loose  and  dense  backfill  materials  used 
in  the  HONDO  calculations  are  all  nonlinear  and  inelastic.  For  the 
SAP  IV  modal  analyses,  linear  elastic  property  specifications  were  re- 
quired. 

Elastic  properties  for  the  lower  bound  case  were  determined  by  cal- 
culating E and  v for  the  free-field  and  loose  backfill  materials 
based  on  (l)  the  loading  secant  constrained  modulus  to  300  psi  (20.7  bars) 
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A. 4 DISCUSSION  OF  RESULTS 

The  in  vacuo  analyses  indicated  that  the  frequency  responses  of  the 
structure  based  on  the  use  of  constant  strain  elements  were  not  appreciably 
different  from  those  obtained  using  higher  order  elements.  They  also  indi- 
cated that  use  of  four-element  thicknesses  across  all  sections  would  pro- 
vide an  adequate  idealization  of  the  structure  for  the  structure /backfill 
interaction  calculations.  The  in  vacuo  results  for  constant  strain  elements 
and  four-element  thicknesses  gave  frequencies  for  the  first  10  structural 
or  deformation  modes  ranging  from  15.08  to  121. 7 Hz. 

The  embedded  analyses  indicated  a significant  effect  of  material  prop- 
erty variation  on  the  frequencies  of  the  buried  structure,  i.e.,  frequencies 
obtained  for  the  first  25  modes  using  lower  bound  properties  ranged  from 
4.315  to  14.35  Hz,  while  those  for  upper  bound  properties  ranged  from  11.99 
to  54.87  Hz.  Embedment  in  either  set  of  soils  obviously  lowers  the  fre- 
quencies from  those  obtained  in  vacuo,  but  it  is  difficult  to  make  mode-for- 
mode  comparisons  because  the  deflected  shapes  are  quite  different,  e.g.. 


from  the  UX  vertical  stress  versus  vertical  strain  relations,  and  (2) 
the  associated  secant  to  the  UX  loading  paths  of  principal  stress  diff- 
erence versus  mean  normal  stress.  The  values  obtained  (along  with  their* 
implied  shear  modulus  G values)  are  given  in  Table  A. 3.  The  first 
25  natural  frequencies  calculated  for  the  structure  embedded  in  materials 
having  lower  bound  elastic  properties  are  listed  in  Table  A. 4.  The  struc- 
ture mode  shapes  associated  with  each  of  the  25  frequencies  are  shown  in 
Figure  A. 5;  mode  shapes  for  both  the  structure  and  the  surrounding  soil  are 
shown  for  the  first  12  of  these  frequencies  in  Figure  A. 6. 

Elastic  properties  for  the  upper  bound  case  were  determined  by  calcu- 
lating E and  v for  the  free-field  and  dense  backfill  materials  based 
on  secant-  to  300-psi  (20. 7- bar)  moduli  from  UX  unloading  relations.  The 
values  obtained  are  given  in  Table  A. 3.  The  first  25  natural  frequencies 
calculated  for  the  upper  bound  case  are  listed  in  Thble  A. 4;  the  associated 
mode  shapes  for  the  structure  are  shown  in  Figure  A. 7. 


compare  Figure  A. 3 with  Figures  A. 5,  A. 6,  and  A. 7.  The  embedded  eigenvalue 
analyses  do  indicate,  however,  that  the  structure /backfill  interaction 
calculations  have  a sufficient  range  of  frequency  response  to  permit  an 
adequate  assessment  of  the  structure  stresses  and  motions  under  the 
postulated  airblast  environments.  The  lowest  frequencies  that  can  be 
fully  transmitted  in  these  calculations  will  be  about  5 Hz;  the  lowest 
modal  frequencies  calculated  from  the  embedded  analyses  were  4.315  Hz 
with  the  lower  bound  properties  and  11.99  Hz  with  the  upper  bound  proper- 
ties. The  upper  cutoff  frequencies  for  the  planned  calculations  are 
assumed  to  be  25  Hz  for  the  loose  backfill  case  and  50  Hz  for  the  case 
of  dense  backfill;  the  embedded  modal  analyses  using  lower  bound  moduli 
gave  a 25th  mode  frequency  of  Just  l4 . 35  Hz,  and  using  upper  bound  moduli, 
the  frequency  for  the  20th  mode  was  49.11  Hz. 
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Able  A.  4 Calculated  Natural  Frequencies  for  Embedded  Analyses. 


Calculated  Calculated 


Mode  Using  Using 

No.  Lower  Bound  Material  Properties  Upper  Bound  Material  Properties 


1 

4.315 

11.99 

2 

4.350 

16.68 

3 

5.054 

17.17 

4 

5.978 

17.43 

5 

6.044 

25.28 

6 

6.825 

25.47 

7 

7.649 

27.80 

8 

8.98  9 

31.07 

9 

9.032 

31.88 

10 

9.189 

34.18 

11 

9.369 

37.92 

12 

10.27 

39.69 

13 

10.52 

40.49 

14 

10.65 

40.95 

15 

10.87 

41.96 

16 

11.09 

42.39 

17 

11.28 

45.80 

18 

11.37 

47.55 

19 

11.87 

47.96 

20 

12.10 

49.11 

21 

12.35 

51.22 

22 

12.64 

51.55 

23 

12.83 

52.50 

24 

13.39 

52.61 

25 

14.35 

54.87 

,n  vacuo  condition 


icies  associated  mode  s napes  from  the  in  /acuo  eigenval 
!8s  structural  approximation  using  constant  strain  elements 


HO  FT 
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Frequencies  and  associated  mode  shapes  of  structure  from  the  embedded  eigenvalue 
calculation  vising  upper  bound  properties. 


APPENDIX  B 

COMPARISON  OF  CALCULATIONAL  RESULTS  WITH 
RIGID  BODY  MOTION  ANALYSIS 

"Hie  objective  of  this  appendix  is  to  compare  the  results  from  a 
rigid  body  vertical  motion  prediction  method  for  shallow-buried  box 
structures ^ with  rigid  body  motion  from  a HO NDO  code  calculation  of  the 
same  problem.  As  a basis  of  comparison,  the  l-MT/300-psi/dense  backfill 
calculation  was  chosen.  The  geometry  of  this  problem  is  shown  in  Fig- 
ure 2.1,  and  the  airblast  loading  is  shown  in  Figure  2.16.  The  methods 
of  Reference  5»  which  are  basically  those  of  a SDF  analysis,  were  used 
to  predict  the  rigid  body  motion  for  a case  that  is  as  close  to  the 
HONDO  code  calculation  as  was  possible  within  the  limitations  of  the 
model. 

B.l  SINGLE-DEGREE-OF-FREEDOM  ANALYSES 

The  problem  analyzed  is  shown  in  Figure  B.l.  The  structure  was 
assumed  to  be  rigid  and  square  in  plan  view.  The  300-psi  airblast 
loading  was  used  for  the  1-MT  weapon  over  the  center  of  the  structure 
(Figure  2.16)  with  a zero  rather  than  a 10-ms  rise  time. 

The  ID  wave-propagation-analysis  hand  calculation  procedures  were 
used  to  compute  the  free-field  vertical  stress  time  histories  at  the 
structure  roof  (1. 52-metre  depth)  and  floor  (9.l4-metre  depth)  levels. 
Figure  B.2  shows  the  free-field  vertical  stress  time  histories  for  the 
roof  and  floor  level. 

The  SDF  equation  of  motion,  neglecting  shear  on  the  sides  of  the 
buried  structure,  is 

Mv3  = A[oR(t)  - op(t)]  (B.l) 


where 

M * mass  of  structure 

v * acceleration  of  structure 
s 
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A = vertically  projected  area  of  structure 
aR(t)  = see  Equation  B.2 
Op(t)  = see  Equation  B.3 

In  Reference  10,  the  vertical  stress  oR  acting  on  the  roof  of  a rec- 
tangular rigid  structure  buried  in  a half-space  of  elastic  soil  is  ap- 
proximated by 


°R  = 2°ffr. 


(B.2) 


where 

= free-field  vertical  stress  at  roof  level 
R 

P1  = mass  density  of  soil  at  roof  level 
CT  = wave  velocity  of  soil  at  roof  level 

h 

vg  = velocity  of  structure 

In  this  same  reference,  the  stress  acting  on  the  floor  is  approximated  by 


v 

s 


(B.3) 


where 

Of.*  = free-field  stress  at  floor  level 
F 

P2  = mass  density  of  soil  at  floor  level 

C = wave  velocity  of  soil  at  floor  level 
2 

The  result  of  combining  Equations  B.l,  B.2,  and  B.3  gives 


v 

s 


M ^2affR  " ° 

n 


Av 


ffp^  M W 


(B.U) 


The  stress  time  history  for  the  net  or  resultant  stress  (2off  - off.  ) is 
given  in  Figure  B.3.  Equation  B.4  was  solved  numerically  using  the  values 
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mm 


of  A 


M 


given  in  Figure  B.7,  and  the  stress  time 


history  in  Figure  B.3  to  determine  the  velocity  time  history  of  the 
structure.  A time  increment  of  2 ms  was  used.  Hie  displacement  time 
history  of  the  structure  was  determined  by  integrating  the  velocity 
time  history. 


B.2  COMPARISONS  WITH  2D  CALCULATION 

RESULTS 

Hie  free- field  stress  and  structure  velocity  and  displacement 
time  histories  from  the  SDF  analyses  are  compared  with  those  from 
the  2D  calculation  for  the  l-MT/300-psi/dense  backfill  calculations 
in  Figures  B.4  through  B,7. 

As  shown  in  Figure  B.4,  the  free-field  vertical  stress  at  the 
1.52-metre  depth  (structure  roof  level)  from  the  ID  analysis  are  com- 
pared with  free-field  vertical  stresses  from  elements  179  and  180, 
which  are  located  at  depths  of  1.0  and  1.9  metres,  respectively,  at  a 
distance  of  2.67  metres  upstream  from  the  structure.  All  three  time 
histories  were  adjusted  to  be  compatible  with  the  arrival  time  of  sur- 
face stress  of  27.5  ms  above  the  center  of  the  structure.  Also,  in 
Figure  B.4,  the  free-field  vertical  stress  time  history  from  the  ID 
analyses  compares  very  well  with  those  from  the  2D  calculations  with 
the  major  discrepancy  being  the  rise  time  to  peak  stress.  The  surface 
airblast  relation  for  the  ID  analyses  contained  a shock  front,  i.e., 
0-ms  rise  time,  while  that  for  the  2D  calculation  contained  a 10-ms 
rise  time. 

In  Figure  B.5,  the  free-field  vertical  stress  at  the  9.l4-metre 
depth  (structure  floor  level)  from  the  ID  analysis  is  compared  with  the 
free-field  vertical  stresses  from  elements  189  and  190,  which  are  lo- 
cated at  depths  of  8.75  and  9.52  metres,  respectively,  at  a distance  of 
2.67  metres  upstream  from  the  structure.  Hie  arrival  times  for  all 
cases  were  also  adjusted  to  correspond  with  an  airblast  time  of  arrival 
at  27.5  ms  at  the  ground  surface  over  the  center  of  the  structure.  As 
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